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The Physical Implementation of Quantum Computation 

David P. DiVincenzo 

IBM T. J. Watson Research Center, Yorktown Heights, NY 10598 USA 
Abstract 



After a brief introduction to the principles and promise of quantum information processing, the require- 
ments for the physical implementation of quantum computation are discussed. These five requirements, 
plus two relating to the communication of quantum information, are extensively explored and related to 
the many schemes in atomic physics, quantum optics, nuclear and electron magnetic resonance spectro- 
2 ° ' scopy, superconducting electronics, and quantum-dot physics, for achieving quantum computing. 



305 i. Introduction 



The advent of quantum information processing, as an abstract concept, has given birth to a 
325 great deal of new thinking, of a very concrete form, about how to create physical comput- 

ing devices that operate in the hitherto unexplored quantum mechanical regime. The efforts 
now underway to produce working laboratory devices that perform this profoundly new 
form of information processing are the subject of this book. 

In this chapter I provide an overview of the common objectives of the investigations 
reported in the remainder of this special issue. The scope of the approaches, proposed and 
underway, to the implementation of quantum hardware is remarkable, emerging from spe- 
;ractions . 355 cialties in atomic physics [1], in quantum optics [2], in nuclear [3] and electron [4] mag- 

netic resonance spectroscopy, in superconducting device physics [5], in electron physics [6], 
and in mesoscopic and quantum dot research [7]. This amazing variety of approaches has 

363 arisen because, as we will see, the principles of quantum computing are posed using the 

most fundamental ideas of quantum mechanics, ones whose embodiment can be contem- 
plated in virtually every branch of quantum physics. 

The interdisciplinary spirit which has been fostered as a result is one of the most pleas- 
ant and remarkable features of this field. The excitement and freshness that has been pro- 
duced bodes well for the prospect for discovery, invention, and innovation in this endeavor. 



2. Why Quantum Information Processing? 

The shortest of answers to this question would be, why not? The manipulation and trans- 
mission of information is today carried out by physical machines (computers, routers, scan- 
ners, etc.), in which the embodiment and transformations of this information can be de- 
scribed using the language of classical mechanics. But the final physical theory of the 
world is not Newtonian mechanics, and there is no reason to suppose that machines follow- 
ing the laws of quantum mechanics should have the same computational power as classical 
machines; indeed, since Newtonian mechanics emerges as a special limit of quantum me- 
chanics, quantum machines can only have greater computational power than classical ones. 
The great pioneers and visionaries who pointed the way towards quantum computers, 
Deutsch [8], Feynman [9], and others, were stimulated by such thoughts. Of course, by a 
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similar line of reasoning, it may well be asked whether machines embodying the principles 
of other refined descriptions of nature (perhaps general relativity or string theory) may have 
even more information processing capabilities; speculations exist about these more exotic 
possibilities, but they are beyond the scope of the present discussion. 

But computing with quantum mechanics really deserves a lot more attention than worm- 
hole computing or quantum-gravity computing; quantum computing, while far in the future 
from the perspective of CMOS roadmaps and projections of chip fab advances, can cer- 
tainly be seen as a real prospect from the perspective of research studies in quantum phy- 
sics. It does not require science fiction to envision a quantum computer; the proposals 
discussed later in this issue paint a rather definite picture of what a real quantum computer 
will look like. 

So, how much is gained by computing with quantum physics over computing with 
classical physics? We do not seem to be near to a final answer to this question, which is 
natural since even the ultimate computing power of classical machines remains unknown. 
But the answer as we know it today has an unexpected structure; it is not that quantum 
tools simply speed up all information processing tasks by a uniform amount. By a stand- 
ard complexity measure (i.e., the way in which the number of computational steps re- 
quired to complete a task grows with the "size" n of the task), some tasks are not sped 
up at all [10] by using quantum tools (e.g., obtaining the nth iterate of a function 
/(/"(• • •/(■*) ■ • •)) H some are sped up moderately (locating an entry in a database of n 
entries [12]), and some are apparently sped up exponentially (Shor's algorithm for factor- 
ing an n-digit number [13]). 

In other types of information processing tasks, particularly those involving communica- 
tion [14], both quantitative and qualitative improvements are seen [15]: for certain tasks 
(choosing a free day for an appointment between two parties from out of n days) there is a 
quadratic reduction of the amount of communicated data required, if quantum states rather 
than classical states are transmitted [16]. For some tasks (the "set disjointness problem", 
related to allocating non-overlapping segments of a shared memory in a distributed compu- 
tation) the reduction of required communication is exponential [17]. Finally, there are tasks 
that are doable in the quantum world that have no counterpart classically: quantum crypto- 
graphy provides an absolute secrecy of communication between parties that is impossible 
classically [18]. And for some games, winning strategies become possible with the use of 
quantum resources that are not available otherwise [19, 20]. 

This issue, and this chapter, are primarily concerned with the "hows" of quantum com- 
puting rather than the "whys," so we will leave behind the computer science after this 
extremely brief mention. There is no shortage of other places to obtain more information 
about these things; I recommend the recent articles by Aharonov [21] and by Cleve [22]; 
other general introductions [23] will give the reader pointers to the already vast specialized 
literature on this subject. 



3. Realizing Quantum Computation 

Let me proceed with the main topic: the physical realization of quantum information pro- 
cessing. As a guide to the remainder of the special issue, and as a means of reviewing the 
basic steps required to make quantum computation work, I can think of no better plan than 
to review a set of basic criteria that my coworkers and I have been discussing over the last 
few years [24] for the realization of quantum computation (and communication), and to 
discuss the application of these criteria to the multitude of physical implementations that 
are found below. 

So, without further ado, here are the 
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Five (plus two) requirements for the implementation of quantum c mputati n 

/. A scalable physical system with well characterized qubits 

For a start, a physical system containing a collection of qubits is needed. A qubit (or, 
more precisely, the embodiment of a qubit) is [25] simply a quantum two-level system like 
the two spin states of a spin 1/2 particle, like the ground and excited states of an atom, or 
like the vertical and horizontal polarization of a single photon. The generic notation for a 
qubit state denotes one state as |0) and the other as |1). The essential feature that distin- 
guishes a qubit from a bit is that, according to the laws of quantum mechanics, the per- 
mitted states of a single qubit fills up a two-dimensional complex vector space; the general 
state is written a |0) + b | 1), where a and b are complex numbers, and a normalization 
convention |a| 2 + |q — 1 is normally adopted. The general state of two qubits, 
a |00) + b |01) + c 1 10) + d |1 1), is a four-dimensional vector, one dimension for each dis- 
tinguishable state of the two systems. These states are generically entangled, meaning that 
they cannot be written as a product of the states of two individual qubits. The general state 
of n qubits is specified by a 2 n -dimensional complex vector. 

A qubit being "well characterized" means several different things. Its physical param- 
eters should be accurately known, including the internal Hamiltonian of the qubit (which 
determines the energy eigenstates of the qubit, which are often, although not always, taken 
as the |0) and |1) states), the presence of and couplings to other states of the qubit, the 
interactions with other qubits, and the couplings to external fields that might be used to 
manipulate the state of the qubit. If the qubit has third, fourth, etc., levels, the computer's 
control apparatus should be designed so that the probability of the system ever going into 
these states is small. The smallness of this and other parameters will be determined by the 
capabilities of quantum error correction, which will be discussed under requirement 3. 

Recognizing a qubit can be trickier than one might think. For example, we might consid- 
er a pair of one-electron quantum dots that share a single electron between them as a two- 
qubit system. It is certainly true that we can denote the presence or absence of an electron 
on each dot by |0) and |1), and it is well known experimentally how to put this system into 
the "entangled" state 1/a/2 (|01) + [ 1 0) ) in which the electron is in a superposition of 
being on the left dot and the right dot. But it is fallacious to consider this as a two-qubit 
system; while the states |00) and |11) are other allowed physical states of the dots, super- 
selection principles forbid the creation of entangled states involving different particle num- 
bers such as l/y/2 (|00) + 

It is therefore false to consider this as a two-qubit system, and, since there are not two 
qubits, it is nonsense to say that there is entanglement in this system. It is correct to say 
that the electron is in a superposition of different quantum states living on the two different 
dots. It is also perfectly correct to consider this system to be the embodiment of a single 
qubit, spanned by the states (in the misleading notation above) |01) ("electron on the right 
dot") and |10) ("electron on the left dot"). Indeed, several of the viable proposals, includ- 
ing the ones by Schon, Averin, and Tanamoto in this special issue, use exactly this system 
as a qubit. However, false lines of reasoning like the one outlined here have sunk various 
proposals before they were properly launched (no such abortive proposals are represented 
in this book, but they can be found occasionally in the literature). 

An amazing variety of realizations of the qubit are represented in this volume. There is a 
very well developed line of work that began with the proposal of Cirac and Zoller [1] for 
an ion-trap quantum computer, in which, in its quiescent state, the computer holds the 
qubits in pairs of energy levels of ions held in a linear electromagnetic trap. Various pairs 
of energy levels (e.g., Zeeman-degenerate ground states, as are also used in the NMR 
approach [3] discussed by Cory) have been proposed and investigated experimentally. The 
many neutral-atom proposals (see chapters by Kimble [2], Deutsch [26], and Briegel 
[27]) use similar atomic energy levels of neutral species. These atomic-physics based propos- 
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als use other auxiliary qubits such as the position of atoms in a trap or lattice, the presence 
or absence of a photon in an optical cavity, or the vibrational quanta of trapped electrons, 
ions or atoms (in the Platzman proposal below [6] this is the primary qubit). Many of the 
solid-state proposals exploit the fact that impurities or quantum dots have well characterized 
discrete energy level spectra; these include the spin states of quantum dots (see chapters by 
Loss [7] and Imamoglu [2]), the spin states of donor impurities (see Kane [4]), and the 
orbital or charge states of quantum dots (see Tanamoto [7]). Finally, there are a variety of 
interesting proposals which use the quantized states of superconducting devices, either ones 
involving the (Cooper-pair) charge (see SchOn, Averin), or the flux (see Moou) [5]. 

2. The ability to initialize the state of the qubits to a simple fiducial state, such as |000 . . .} 
This arises first from the straightforward computing requirement that registers should be 
initialized to a known value before the start of computation. There is a second reason for 
this initialization requirement: quantum error correction (see requirement 3 below) requires a 
continuous, fresh supply of qubits in a low-entropy state (like the |0) state). The need for a 
continuous supply of Os, rather than just an initial supply, is a real headache for many pro- 
posed implementations. But since it is likely that a demonstration of a substantial degree of 
quantum error correction is still quite some time off, the problem of continuous initialization 
does not have to be solved very soon; still, experimentalists should be aware that the speed 
with which a qubit can be zeroed will eventually be a very important issue. If the time it 
takes to do this initialization is relatively long compared with gate-operation times (see 
requirement 4), then the quantum computer will have to equipped with some kind of "qubit 
conveyor belt", on which qubits in need of initialization are carried away from the region in 
which active computation is taking place, initialized while on the "belt", then brought back 
to the active place after the initialization is finished. (A similar parade of qubits will be 
envisioned in requirement 5 for the case of low quantum-efficiency measurements [28].) 

There are two main approaches to setting qubits to a standard state: the system can either 
be "naturally" cooled when the ground state of its Hamiltonian is the state of interest, or 
the standard state can be achieved by a measurement which projects the system either into 
the state desired or another state which can be rotated into it. These approaches are not 
fundamentally different from one another, since the projection procedure is a form of cool- 
ing; for instance, the laser cooling techniques used routinely now for the cooling of ion 
states to near their ground state in a trap [1] are closely connected to the fluorescence 
techniques used to measure the state of these ions. A more "natural" kind of cooling is 
advocated in many of the electron spin resonance based techniques (using quantum dots or 
impurities) [7, 4] in which the spins are placed in a strong magnetic field and allowed to 
align with it via interaction with their heat bath. In this kind of approach the time scale will 
be a problem. Since the natural thermalization times are never shorter than the decoherence 
time of the system, this procedure will be too slow for the needs of error correction and a 
"conveyor belt" scheme would be required. Cooling by projection, in which the Hamilto- 
nian of the system and its environment are necessarily perturbed strongly, will have a time 
scale dependent on the details of the setup, but potentially much shorter than the natural 
relaxation times. One cannot say too much more at this point, as the schemes for measure- 
ment have in many cases not been fully implemented (see requirement 5). In the NMR 
quantum computer implementations to date (see Cory below), cooling of the initial state 
has been foregone altogether; it is acknowledged [3] that until some of the proposed cool- 
ing schemes are implemented (a nontrivial thing to do), NMR can never be a scalable 
scheme for quantum computing. 

3. Long relevant decoherence times, much longer than the gate operation time 
Decoherence times characterize the dynamics of a qubit (or any quantum system) in contact 
with its environment. The (somewhat overly) simplified definition of this time is that it is 
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the characteristic time for a generic qubit state \y/) = a |0) + b |1) to be transformed into 
the mixture p = \a\ 2 |0) (0| + \b\ 2 |1) (1|. A more proper characterization of decoherence, in 
which the decay can depend on the form of the initial state, in which the state amplitudes 
may change as well, and in which other quantum states of the qubit can play a role (in a 
special form of state decay called "leakage" in quantum computing [29, 30]), is rather more 
technical than I want to get here; but see Refs. [31] and [32] for a good general discussion 
of all these. Even the simplest discussion of decoherence that I have given here should also 
be extended to include the possibility that the decoherence of neighboring qubits is corre- 
lated. It seems safest to assume that they will be neither completely correlated nor comple- 
tely uncorrelated, and the thinking about error correction has taken this into account. 

Decoherence is very important for the fundamentals of quantum physics, as it is identi- 
fied as the principal mechanism for the emergence of classical behavior. For the same rea- 
son, decoherence is very dangerous for quantum computing, since if it acts for very long, 
the capability of the quantum computer will not be so different from that of a classical 
machine. The decoherence time must be long enough that the uniquely quantum features of 
this style of computation have a chance to come into play. How long is "long enough" is 
also indicated by the results of quantum error correction, which I will summarize shortly. 

I have indicated that the "relevant" decoherence times should be long enough. This em- 
phasizes that a quantum particle can have many decoherence times pertaining to different 
degrees of freedom of that particle. But many of these can be irrelevant to the functioning 
of this particle as a qubit. For example, the rapid decoherence of an electron's position state 
in a solid state environment does not preclude its having a very long spin coherence time, 
and it can be arranged that this is the only time relevant for quantum computation. Which 
time is relevant is determined by the choice of the qubit basis states |0) and |1); for exam- 
ple, if these two states correspond to different spin states but identical orbital states, then 
orbital decoherence will be irrelevant. 

One might worry that the decoherence time necessary to do a successful quantum com- 
putation will scale with the duration of the computation. This would place incredibly strin- 
gent requirements on the physical system implementing the computation. Fortunately, in 
one of the great discoveries of quantum information theory (in 1995—1996), it was found 
that error correction of quantum states is possible [33] and that this correction procedure 
can be successfully applied in quantum computation [34], putting much more reasonable 
(although still daunting) requirements on the needed decoherence times. 

In brief, quantum error correction starts with coding; as in binary error correction codes, 
in which only a subset of all boolean strings are "legal" states, quantum error correction 
codes consist of legal states confined to a subspace of the vector space of a collection of 
qubits. Departure from this subspace is caused by decoherence. Codes can be chosen such 
that, with a suitable sequence of quantum computations and measurements of some ancil- 
lary qubits, the error caused by decoherence can be detected and corrected. As noted above, 
these ancillary qubits have to be continuously refreshed for use. I will not go much farther 
into the subject here, see [31] for more. It is known that quantum error correction can be 
made fully fault tolerant, meaning that error correction operations can be successfully inter- 
mingled with quantum computation operations, that errors occurring during the act of error 
correction, if they occur at a sufficiently small rate, do no harm, and that the act of quan- 
tum computation does not itself cause an unreasonable proliferation of errors. 

These detailed analyses have indicated the magnitude of decoherence time scales that are 
acceptable for fault-tolerant quantum computation. The result is that, if the decoherence 
time is 10 4 — 10 5 times the "clock time" of the quantum computer, that is, the time for the 
execution of an individual quantum gate (see requirement 4), then error correction can be 
successful. This is, to tell the truth, a rather stringent condition, quantum systems frequently 
do not have such long decoherence times. But sometimes they do, and .our search for a 
successful physical implementation must turn towards these. At least this result says that 
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the required decoherence rate does not become ever smaller as the size and duration of the 
quantum computation grows. So, once the desired threshold is attainable, decoherence will 
not be an obstacle to scalable quantum computation. 

Having said this, it must be admitted that it will be some time before it is even possible 
to subject quantum error correction to a reasonable test. Nearly all parts of requirements 
1-5 must be in place before such a test is possible. And even the most limited application 
of quantum error correction has quite a large overhead: roughly 10 ancillary qubits must be 
added for each individual qubit of the computation. Fortunately, this overhead ratio grows 
only logarithmically as the the size of the quantum computation is increased. 

In the short run, it is at least possible to design and perform experiments which measure 
the decoherence times and other relevant properties (such as the correlation of decoherence 
of neighboring qubits) of candidate implementations of qubits. With such initial test experi- 
ments, caution must be exercised in interpreting the results, because decoherence is a very 
system-specific phenomenon, depending on the details of all the qubits' couplings to var- 
ious environmental degrees of freedom. For example, the decoherence time of the spin of 
an impurity in the bulk of a perfect semiconductor may not be the same as its decoherence 
time when it is near the surface of the solid, in the immediate neighborhood of device 
structures designed to manipulate its quantum state. Test experiments should probe decoher- 
ence in as realistic a structure as is possible. 

4. A "universal" set of quantum gates 
This requirement is of course at the heart of quantum computing. A quantum algorithm is 
typically specified [8] as a sequence of unitary transformations £/,, U 2 , U 3 , . . each acting 
on a small number of qubits, typically no more than three. The most straightforward tran- 
scription of this into a physical specification is to identify Hamiltonians which generate 
these unitary transformations, viz., £/, = e iHlt ^ 9 U 2 = <? ,7, 2'/ fi , £/ 3 = e iH ^^ h , etc.; then, the 
physical apparatus should be designed so that //, can be turned on from time 0 to time t, 
then turned off and H 2 turned on from time t to time It, etc. 

Would that life were so simple! In reality what can be done is much less, but much less 
can be sufficient. Understanding exactly how much less is still enough, is the main compli- 
cation of this requirement. In all the physical implementations discussed in this volume, 
only particular sorts of Hamiltonians can be turned on and off; in most cases, for example! 
only two-body (two-qubit) interactions are considered. This immediately poses a problem 
for a quantum computation specified with three-qubit unitary transformations; fortunately, 
of course, these can always be re-expressed in terms of sequences of one- and two-body 
interactions [35], and the two-body interactions can be of just one type [36], the "quantum 
XOR" or "cNOT". There are some implementations in which multi-qubit gates can be 
implemented directly [37]. 

However, this still leaves a lot of work to do. In some systems, notably in NMR (see 
Cory), there are two-body interactions present which cannot be turned off, as well as others 
which are switchable. This would in general be fatal for quantum computation, but the 
particular form of the fixed interactions permit their effects to be annulled by particular 
"refocusing" sequences of the controllable interactions, and it has recently been discovered 
[38] that these refocusing sequences can be designed and implemented efficiently. 

For many other systems, the two-body Hamiltonian needed to generate directly the cNOT 
unitary transformation is not available. For example, in the quantum-dot proposal described 
by Loss below [7], the only two-body interaction which should be easily achievable is the 
exchange interaction between neighboring spins, H oc S t - S i+l ; in the Imamoglu chapter [2] 
the attainable interaction is of the XY type, i.e., H oc S^ y , + S iy S jr An important observa- 
tion is that with the appropriate sequence of exchange or XY interactions, in conjunction 
with particular one-body interactions (which are assumed to be more easily doable), the 
cNOT transformation can be synthesized [39]. It is incumbent on 'each implementation pro- 
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posal to exhibit such a sequence for producing the cNOT using the interactions that are 
naturally realizable. 

Often there is also some sophisticated thinking required about the time profile of the 
two-qubit interaction. The naive description above uses a "square pulse" time profile, but 
often this is completely inappropriate; for instance, if the Hamiltonian can also couple the 
qubit to other, higher-lying levels of the quantum system, often the only way to get the 
desired transformation is to turn on and off the interaction smoothly and slowly enough that 
an adiabatic approximation is accurate [29, 30] (in a solid-state context, see also [40]). The 
actual duration of the pulse will have to be sufficiently long that any such adiabatic require- 
ment is satisfied; then typically only the time integral \dtH{t) is relevant for the quantum 
gate action. The overall time scale of the interaction pulse is also controlled by the attain- 
able maximum size of the matrix elements of H(t), which will be determined by various 
fundamental considerations, like the requirement that the system remains in the regime of 
validity of a linear approximation, and practical considerations, like the laser power that can 
be concentrated on a particular ion. Given these various constraints, the "clock time" of the 
quantum computer will be determined by the time interval needed such that two consecu- 
tive pulses have negligible overlap. 

Another consideration, which does not seem to present a problem with any current im- 
plementation schemes, but which may be an issue in the future, is the classicality of the 
control apparatus. We say that the interaction Hamiltonian H(t) has a time profile which is 
controlled externally by some "classical" means, that is, by the intensity of a laser beam, 
the value of a gate voltage, or the current level in a wire. But each of these control devices 
is made up themselves of quantum mechanical parts. When we require that these behave 
classically, it means that their action should proceed without any entanglement developing 
between these control devices and the quantum computer. Estimates indicate that this entan- 
glement can indeed be negligible, but this effect needs to be assessed for each individual 
case. 

In many cases it is impossible to turn on the desired interaction between a pair of qubits; 
for instance, in the ion-trap scheme, no direct interaction is available between the ion-level 
qubits [1]. In this and in other cases, a special quantum subsystem (sometimes referred to 
as a "bus qubit") is used which can interact with each of the qubits in turn and mediate the 
desired interaction: for the ion trap, this is envisioned to be the vibrational state of the ion 
chain in the trap; in other cases it is a cavity photon whose wavefunction overlaps all the 
qubits. Unfortunately, this auxiliary quantum system introduces new channels for the envi- 
ronment to couple to the system and cause decoherence, and indeed the decoherence occur- 
ring during gate operation is of concern in the ion-trap and cavity-quantum electrodynamics 
schemes. 

Some points about requirement 4 are important to note in relation to the implementation 
of error correction. Successful error correction requires fully parallel operation, meaning 
that gate operations involving a finite fraction of all the qubits must be doable simulta- 
neously. This can present a problem with some of the proposals in which the single "bus 
qubit" is needed to mediate each interaction. On the other hand, the constraint that interac- 
tions are only among nearest neighbors in a lattice, as in many of the solid-state proposals, 
does allow for sufficient parallelism [41]. 

Quantum gates cannot be implemented perfectly; we must expect both systematic and 
random errors in the implementation of the necessary Hamiltonians. Both types of errors 
can be viewed as another source of decoherence and thus error correction techniques are 
effective for producing reliable computations from unreliable gates, if the unreliability is 
small enough. The tolerable unreliability due to random errors is in the same vicinity as the 
decoherence threshold, that is, the magnitude of random errors should be 10" 4 — 10" 5 per 
gate operation or so. It might be hoped that systematic errors could be virtually eliminated 
by careful calibration; but this will surely not always be the case. It seems harder to give a 
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good rule for how much systematic error is tolerable, the conservative estimates give a 
very, very small number (the square of the above) [31], but on the other hand there seems 
to be some evidence that certain important quantum computations (e.g., the quantum Four- 
ier transform) can tolerate a very high level of systematic error (over- or under-rotation). 
Some types of very large errors may be tolerable if their presence can be detected and 
accounted for on the fly (we are thinking, for example, about charge switching in semicon- 
ductors or superconductors). 

Error correction requires that gate operations be done on coded qubits, and one might 
worry that such operations would require a new repertoire of elementary gate operations for 
the base-level qubits which make up the code. For the most important error correction 
techniques, using the so called "stabilizer" codes, this is not the case. The base-level toolkit 
is exactly the same as for the unencoded case: one-bit gates and cNOTs, or any gate reper- 
toire that can produce these, are adequate. Sometimes the use of coding can actually reduce 
the gate repertoire required: in the work on decoherence free subspaces and subsystems, 
codes are introduced using blocks of three and four qubits for which two-qubit exchange 
interactions alone are enough to implement general quantum computation [42, 43]. This 
simplification could be very useful in the quantum-dot [7] or semiconductor-impurity [4] 
implementations. 

5. A qubit-specific measurement capability 
Finally, the result of a computation must be read out, and this requires the ability to 
measure specific qubits. In an ideal measurement, if a qubit's density matrix is 
p=p\0) <0| + (1 -p) |1) <1| + a |0> <1| + a* |1> (01, the measurement should give out- 
come "0" with probability p and "1" with probability 1 -p independent of a and of any 
other parameters of the system, including the state of nearby qubits, and without changing 
the state of the rest of the quantum computer. If the measurement is "non-demolition", that 
is, if in addition to reporting outcome "0" the measurement leaves the qubit in state |0) <0|, 
then it can also be used for the state preparation of requirement 2; but requirement 2 can be 
fulfilled in other ways. 

Such an ideal measurement as I have described is said to have 100% quantum efficiency* 
real measurements always have less. While the fidelity of a quantum measurement is not 
captured by a single number, the single quantum-efficiency parameter is often a very useful 
way to summarize it, just as the decoherence time is a useful if incomplete summary of the 
damage caused to a quantum state by the environment. 

While quantum efficiency of 100% is desirable, much less is needed for quantum com- 
putation; there is, in fact, a tradeoff possible between quantum efficiency and other re- 
sources which results in reliable computation. As a simple example, if the quantum effi- 
ciency is 90%, then, in the absence of any other imperfections, a computation with a 
single-bit output (a so-called "decision problem", common in computer science) will have 
90% reliability. If 97% reliability is needed, this can just be achieved by rerunning the 
calculation three times. Much better, actually, is to "copy" the single output qubit to three, 
by applying two cNOT gates involving the output qubit and two other qubits set to |0),' 
and measuring those three. (Of course, qubits cannot be "copied", but their value in a 
particular basis can.) In general, if quantum efficiency q is available, then copying to some- 
what more than \/q qubits and measuring all of these will result in a reliable outcome. So, 
a quantum efficiency of 1% would be usable for quantum computation, at the expense of 
hundreds of copies/remeasures of the same output qubit. (This assumes that the measure- 
ment does not otherwise disturb the quantum computer. If it does, the possibilities are 
much more limited.) 

Even quantum efficiencies much, much lower than 1% can be and are used for success- 
ful quantum computation: this is the "bulk" model of NMR (see Cory and [3]), where 
macroscopic numbers of copies of the same quantum computer (different molecules in solu- 
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tion) run simultaneously, with the final measurement done as an ensemble average over the 
whole sample. These kinds of weak measurements, in which each individual qubit is hardly 
disturbed, are quite common and well understood in condensed-matter physics. 

If a measurement can be completed quickly, on the timescale of 10~ 4 of the decoherence 
time, say, then its repeated application during the course of quantum computation is valu- 
able for simplifying the process of quantum error correction. On the other hand, if this fast 
measurement capability is not available, quantum error correction is still possible, but it 
then requires a greater number of quantum gates to implement. 

Other tradeoffs between the complexity and reliability of quantum measurement vs. those 
of quantum computation have recently been explored. It has been shown that if qubits can 
be initialized into pairs of maximally entangled states, and two-qubit measurements in the 
so-called Bell basis (W ± = jOl) ± ) 10), 0 ± = (00) ± |11» are possible, then no two-qubit 
quantum gates are needed, one-bit gates alone suffice [44]. Now, often this tradeoff will not 
be useful, as in many schemes a Bell measurement would require two-bit quantum gates. 

But the overall message, seen in many of our requirements, is that more and more, the 
theoretical study of quantum computation has offered a great variety of tradeoffs for the 
potential implementations: if X is very hard, it can be substituted with more of Y. Of 
course, in many cases both X and Y are beyond the present experimental state of the art; 
but a thorough knowledge of these tradeoffs should be very useful for devising a rational 
plan for the pursuit of future experiments. 



4. Desiderata for Quantum Communication 

For computation alone, the five requirements above suffice. But the advantages of quantum 
information processing are not manifest solely, or perhaps even principally, for straightfor- 
ward computation only. There are many kinds of information-processing tasks, reviewed 
briefly at the beginning, that involve more than just computation, and for which quantum 
tools provide a unique advantage. 

The tasks we have in mind here all involve not only computation but also communica- 
tion. The list of these tasks that have been considered in the light of quantum capabilities, 
and for which some advantage has been found in using quantum tools, is fairly long and 
diverse: it includes secret key distribution, multiparty function evaluation as in appointment 
scheduling, secret sharing, and game playing [14]. 

When we say communication we mean quantum communication: the transmission of 
intact qubits from place to place. This obviously adds more features that the physical appa- 
ratus must have to carry out this information processing. We formalize these by adding two 
more items to the list of requirements: 

6. The ability to intercom ert stationary and flying qubits 

7. The ability to faithfully transmit flying qubits between specified locations 

These two requirements are obviously closely related, but it is worthwhile to consider them 
separately, because some tasks need one but not the other. For instance, quantum cryptogra- 
phy [18] involves only requirement 7; it is sufficient to create and detect flying qubits 
directly. 

I have used the jargon "flying qubits" [2], which has become current in the discussions 
of quantum communication. Using this term emphasizes that the optimal embodiment of 
qubits that are readily transmitted from place to place is likely to be very different from the 
optimal qubits for reliable local computation. Indeed, almost all proposals assume that 
photon states, with the qubit encoded either in the polarization or in the* spatial wavefunc- 
tion of the photon, will be the flying qubit of choice, and indeed, the well developed tech- 
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rrSiJ llght f tran u s ! TUSSIon trough optical fibers provides a very promising system for the 

SSbST^S*?*- 1 W ° U | d n °u' th ° Ugh ' th3t my C0 " eagUeS and 1 have raised th e 
flying qubit U4 e i5 C ] tr ° nS travehn 8 thou S h solids could P rovide another realization of the 

Only a few completely developed proposals exist which incorporate requirements 6 and 
/. ut course, there are a number of quite detailed studies of 7, in the sense that experi- 
ments on quantum cryptography have been very concerned with the preservation of the 
photon quantum state during transmission through optical fibers or through the atmosphere 
However, these studies are rather disconnected from the other concerns of quantum comput- 
ing. Requirement 6 is the really hard one; to date the only theoretical proposal sufficiently 
concrete that experiments addressing it have been planned is the scheme produced by Kim- 
ble and coworkers [46] for unloading a cavity photon into a traveling mode via atomic 
fhe C lZ C h Py ' a " d 1 loadin gj t b y *e time-reversed process. Other promising concepts, like 
tne launching of electrons from quantum dots into quantum wires such that the spin coher- 
ence or the electrons is preserved, need to be worked out more fully. 



5. Summary 

fn'J? 31 is " winnin 8" technology going to be? I don't think that any living mortal has 
an answer to this question, and at this point it may be counterproductive even to ask it 
5£1 • T haVe J ived with <l uan tuni mechanics for a century, our study of quantum 
effects in complex artificial systems like those we have in mind for quantum computing Z 
t 1*/ °" e Ca " SCe h0W or whether a11 the requirements above can be fulfilled 
su^estedhv ZT trade0ff \ n0t envisioned in our present theoretical discussions but 
suggested by further experiments, that might take our investigations in an entirely new path 
Indeed the above discussion, and the other chapters of this special issue, really do not 
cover all the foreseeable approaches. I will mention two of which I am aware- first an- 
other computational paradigm, that of the cellular automaton, is potentially available for 
exploitation. This ,s distinguished from the above "general purpose" approach i^ that k 

lu ion^ e It eVe T blt thr ° Ug , hOUt ** C ° mpUter Wil1 be -hjectedTo the same evo- 

lution rule It is known that general-purpose computation is performable, although with 
considerable overhead, by a cellular automaton. This is true as well for the quantum ver 
sic-n of me cellular automaton, as Lloyd [47] indicated in his original work. New theoreti- 

nermTL ^ [481 f 0 "" ^ eXpHdtl y how relatively simple local rules would 

permit the implementation of some quantum computations. This could point us perhaps 

SUSS .TV 0 * ° f P ° lymer Wkh 3 Strin S of < ubits on * backbone** can be ad* 
uS^SLSSfl? \ 1 \ S P eCt I? SC0 P ic fashion - Experiments are not oriented towards this at 

.J the™ iure. ^ different ' a " d 1 d ° n,t believe k should be excluded 

Second, even more speculative, but very elegant, is the proposal of Kitaev [49] to use 
ST/Y^ WUh of topological excitations, for example noiabelian 

anyons for quantum computing. It is hard to see at the moment how to turn this exciting 
proposal into an experimental program, as no known physical system is agreed to have the 
appropnate topological excitations. But further research in, for example, L quantum Hill 
effect might reveal such a system; more likely, perhaps, is that further understanding ™f "hi 
approach, and that of Freedman and his colleagues [50], will shed more light on doing 
quantum computing using the "standard" approach being considered in this book 8 
Jr^H 0 "™ fl ° ne tWn ! : thC ideaS ° f quantum information theory will continue to 
STrf^ST ' T CC ° n thC fUIther investi g ati °n of the fundamental quantum proper 
ties of complex quantum systems, and will stimulate many creative and exciting develop- 
ments for many years to come. . ueveiop 
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Anomalous periodicity of the current-phase relationship of grain-boundary Josephson junctions 

in high-7 c superconductors 
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The current-phase relation (CPR) for asymmetric 45° Josephson junctions between two rf-wave supercon- 
ductors has been predicted to exhibit an anomalous periodicity. We have used the single-junction interferom- 
eter to investigate the CPR for these kinds of junctions in YBa 2 Cu 3 07_ x thin films. A remarkable amplitude of 
the 7r-periodical component of the CPR has been experimentally found, providing an additional source of 
evidence for the t/-wave symmetry of the pairing state of the cuprates. [SO 163- 1829(99)05629-5] 



A number of experimental results confirm ^2-^2-wave 
symmetry of the pairing state of high-temperature 
superconductors. 1 An unconventional pairing state requires 
the existence of zeros of the order parameter in certain di- 
rections in momentum space. Thermodynamic and spectro- 
scopic measurements do indeed suggest their existence, but 
by themselves they do not exclude conventional s-wave pair- 
ing with nodes. 1 Direct evidence for the c/-wave pairing state 
is provided by phase-sensitive experiments, which are based 
on the Josephson effect. 2 Quite generally, the current-phase 
relationship (CPR) of a Josephson junction, I(<p) is an odd 
periodic function of <p with a period 27r. 3 . Therefore 7(<p) 
can be expanded in a Fourier series 

I{<p) — I\ sin <p + f 2 sin 2 9 H . (1) 

In the tunnel limit we can restrict ourselves to the first two 
terms in Eq. (1). Since the order parameter is bound to the 
crystal lattice, I(<p) of a weak link depends on the orientation 
of the J-wave electrodes with respect to their boundary. The 
existing phase-sensitive experiments exploit possible sign 
changes of l { between different geometries. 2 In this work we 
present a phase-sensitive experimental test of the pairing 
state symmetry of cuprates. Namely, in certain geometries, 
the l x term should vanish by symmetry. In such cases, the 
CPR should exhibit an anomalous periodicity. 

Let us analyze the angular dependence of I\ 2 m a junction 
between two macroscopically tetragonal c/-wave supercon- 
ductors. As emphasized in Ref. 4, also heavily twinned 
orthorhombic materials such as YBa 2 Cu 3 0 7 _ x belong to this 
class, if the twin boundaries have odd symmetry. We con- 
sider an ideally flat interface between two superconducting 
electrodes. Let 0 { (0 2 ) denote the angle between the normal 
to the grain boundary and the a axis in electrode 1 (2), see 
Fig. 1 . If we only keep the lowest-order angular harmonics, 
the symmetry of the problem dictates that 4 

/ 1 = / c cos2^ l cos 2 B 2 + I s sin 2 6 X sin2 0 2 . (2) 

The coefficients I c J 5 are functions of the barrier strength, 
temperature T, etc. The I 2 term results from higher-order 
tunneling processes and we neglect its weak angular depen- 



dence. It is seen from Eq. (2) that the criterion for the obser- 
vation of an anomalous period of the CPR, I x — 0, is realized 
for an asymmetric 45° junction, i.e., a junction with d x 
= 45° and 0 2 = 0. 

The I 2 term is also present in weak links based on con- 
ventional s-wave superconductors but for all known types of 
weak links \I 2 /I\\<\. For instance, for a tunnel junction 
\I 2 il x \< 1 . For a superconductor-normal-metal- 
superconductor (SNS) junction, I<*$\n<pl2 at 7*=0, 5 and the 
Fourier expansion of Eq, (1) leads to I 2 II X --2IS. There- 
fore a possible experimental observation of \I 2 II X \> 1 in an 
asymmetric 45° junction provides direct evidence of d-wave 
symmetry of the pairing state in the cuprates. 

We have investigated the CPR of YBa 2 Cu307-^ thin-film 
bicrystals with asymmetric 45° [001]-tilt grain boundaries as 
sketched in Fig. 1, using a single-junction interferometer 
configuration in which the Josephson junction is inserted into 
a superconducting loop with a small inductance I. In a sta- 
tionary state without fluctuations, the phase difference <p 
across the junction is controlled by applying an external 
magnetic flux 4> e penetrating the loop: <p= <p e - Pf(cp). Here 
<p e = 27r<t> e /<& 0 ; <£> 0 =2.07X 10" 15 Tm 2 is the flux quantum; 
f(<p) = I(<p)/Io is the CPR normalized to the maximal Jo- 
sephson current 7 0 , and f3~27rLI 0 /<£> 0 is the normalized 
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FIG. 1 . Washer-shaped interferometer with one short Josephson 
junction (not in scale). Dimensions are given in the text. 
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critical current. In order to obtain the CPR for the complete 
phase range -tt^^^tt the condition f3<\ has to be ful- 
filled, because for fi> 1 the curve <pO e ) becomes multival- 
ued. Following Ref. 3, we express the effective inductance of 
the interferometer using the derivative /' with respect to <p 
as L int = L[\ + \/pf((p)]. The inductance can be probed by 
coupling the interferometer to a tank circuit with inductance 
L Tt quality factor Q, and resonance frequency a> 0 through 
the mutual inductance M. 8 External flux in the interferometer 
is produced by a current I dc +I r f in the tank coil and can be 
expressed as (p e =27r(! dc + I rf )M/<t> 0 =<p dc + <p rf , where 
M 2 ~k 2 LL T with k a coupling coefficient. Taking into ac- 
count the quasiparticle current in the presence of a voltage V 
across the junction the phase difference is given by the rela- 
tion (p^<p dc +q> rf -pf((p)-27rr((p)V/^ 0 , where r(<p) 
= L/R j(<p) with R j((p) the resistance of the junction. In the 
small-signal limit <p r j<\ and in the adiabatic case 6>r<^l, 
keeping only the first-order terms, the effective inductance 
L ef f of the tank circuit-interferometer system is 



Lint) \ 1 



Thus the phase angle a between the driving current and the 
tank voltage U at the resonance frequency of the tank circuit 
(o 0 is 



tan a(<p)~ 



k 2 QPf'(<p) 



(3) 



Using the relation [1 + Pf f (<p)]d<p = d<p dc which is valid for 
(p r f<\ and q)t<Z\, one can find the CPR from Eq. (3) by 
numerical integration. 

The advantage of the CPR measurement of an asymmetric 
45° junction with respect to the by-now standard phase- 
sensitive tests of pairing symmetry based on the angular de- 
pendence of I } is twofold. First, it avoids the complications 
of the analysis of experiments caused by the presence of the 
term I s . 4 Second, flux trapped in the interferometer washer 
(see Fig. 1) does not invalidate the conclusions about the 
ratio \I 2 fI\ \ and hence about the pairing symmetry, which is 
not the case in standard phase-sensitive tests of the */-wave 
symmetry of the pairing state. 9 

The films of 100-nm thickness were fabricated using stan- 
dard pulsed laser deposition on (001) oriented SrTi0 3 bi- 
crystalline substrates with asymmetric [001] tilt misorienta- 
tion angles of 45° ±1°. The films were subsequently pat- 
terned by Ar ion-beam etching into 4X4-mm 2 square 
washer single-junction interferometer structures (Fig. 1). The 
widths of the junctions were 1-2 /an. The square washer 
holes had a side length of 50 /un. This geometry of the 
interferometer gives L^80 pH. The resistance of a similar 
single junction (without interferometer loop) was measured 
directly and R j> 1 ft was found. Therefore the condition 
for the adiabatic limit <dt< 1 is satisfied. For measurements 
of a(<p dc ) y several tank circuits with inductances 0.2-0.8 /iH 
and resonance frequencies 16-35 MHz have been used. The 
unloaded quality factor of the tank circuits 10<Q< 150 has 
been measured at various temperatures. The coupling factor 
k was determined from the period &I dc of ot{I dc ) using 
MAI dc -<& 0 . Its value varied between 0.03 and 0.09. The 
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FIG. 2. Left panel: Phase angle between the driving current and 
the output voltage measured for sample No. 1 at different tempera- 
tures as a function of the dc current I dc . The curves are shifted 
along the y axis and the data for 7=30 and 40 K are multiplied by 
factor 4 for clarity. From top to bottom, the data correspond to 7 
= 4.2, 10, 15, 20, 30, and 40 K. The data measured for 36° bicrys- 
tals (0^36°, 0 2 ^O) at 7=40 K in the same washer geometry are 
shown for comparison (open circles). Right panel: The same for 
sample No. 3. From top to bottom, the data correspond to 7=4.2, 
10, 15, and 20 K. 

amplitude of I r f was set to produce a flux in the interferom- 
eter smaller than 0.1 <£ 0 t0 ensure the small-signal limit. 

The measurements have been performed in a gas-flow 
cryostat with a five-layer magnetic shielding in the tempera- 
ture range 4.2^7<90 K. The experimental setup was cali- 
brated by measuring interferometers of the same size with 
24° and 36° grain boundaries. We have studied six samples, 
out of which for four samples the 7r-periodic component of 
7(<p) was experimentally observed. At low temperatures for 
two samples (Nos. 1 and 2) the value of I 2 is larger than I x . 
For sample Nos. 3 and 4 I 2 is approximately 10-20% of/! 
and for sample Nos. 5 and 6 7 2 is negligible. As an example 
we plot the phase angle a as a function of the dc current I dc 
for sample Nos. 1 and 3 (Fig. 2). The behavior of sample No. 
1 at low temperatures is defined by the it periodic compo- 
nent of I(<p). The curves for sample No. 3 are 27r-periodic, 
nevertheless for the curve at 7=4.2 K the local minima 
clearly show the presence of a ir-periodic component. 

In order to determine the CPR we assume that the period 
of a(1 dc ) at 7=40 K and A/ rfc = 9.6 /iA, corresponds to 
&(p dc = 2TT. We take <p dc =0 at a maximum or minimum of 
a. This is necessary in order to satisfy l(<p = 0) = 0, as re- 
quired by general principles. 3 The experimentally observed 
shift of the first extreme of a(f dc ) from I dc =0 (Fig. 2) can 
be due to flux trapped in the interferometer washer. Most 
probably, this flux resides in the long junction originated by 
the grain boundary crossing the washer of the interferometer. 
This long junction does not play an active role because the 
Josephson penetration depth is much smaller than the junc- 
tion length, and external fields produced by I dc are smaller 
than the first critical field. Nevertheless, the long junction 
sets the phase difference for I dc ~0 at the small junction. 

In Fig. 3, we show the CPR determined from the data in 
Fig. 2. For all curves we have performed a minimal neces- 
sary shift consistent with 7(<p = 0) = 0. Thus we have as- 
sumed that at <p dc = 0 a minimum of a(<p dc ) is realized. For 
an interferometer with a conventional 5- wave weak link (and 
also for the 36° junction), at ip dc = 0 one gets a maximum of 
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FIG. 3. (a) Josephson current through the junction for sample 
No. 1 as a function of the phase difference <p, determined from the 
data in Fig. 2. The scattering of a(<p) values was reduced by fold- 
ing the data back to the interval (0,7r) and taking the average. Here, 
the symmetry a(<p) = a( — <p) was assumed, (b) The same for 
sample No. 3. The data for the asymmetric 36° bicrystal at 7=40 K 
(open circles) are also shown. 

a{<p dc ). Note that the minimum of a(<p dc ) at <p<j c = 0 implies 
a paramagnetic response of the interferometer in the limit of 
small applied fields. 

The amplitude of the 7r-periodic component of the CPR 
decreases drastically with increasing temperature, and at T 
= 40 K its contribution is negligible for all samples. The 
temperature dependence of I\ and I 2 could be determinated 
with acceptable accuracy for sample No. 1 only. With de- 
creasing T, \I 2 \ grows monotonically down to 7=4.2 K, 
while the I] component exhibits only a weak temperature 
dependence (Fig. 4). 
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FIG. 4. Temperature dependence of the Fourier expansion coef- 
ficients /| 2 determined from the experimental data in Fig. 3(a). 
Solid lines are the Fourier expansion coefficients for the numerical 
data in Fig. 5. 
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FIG. 5. I((p) calculated according to Eq. (64) of Ref. 11 for a 
junction with 0, = 45.5°, 0 2 = O, \d= 1.5, k=0.5, and f c = 60 K. 
I(<p) at 7=40 K for the 36° bicrystal (open circles) was calculated 
with the same parameters except for 0 I = 36°. 

Our experimental results can be understood as follows. It 
is well known that the microstructural properties of the grain 
boundaries, especialy 45° boundaries, are defined by their 
faceted nature. Faceting is an intrinsic property of the grain 
boundaries, 6,7 and, due to rf-wave symmetry of the order pa- 
rameter, the properties of the junctions strongly depend on 
the particular distribution of the facets. Small deviations 
from the ideal geometry of the asymmetric 45° junction lead 
to a finite value of I x . Thus for nearly ideal junctions 
|/ 2 //il>l at T-vO. The region T-T c can be analyzed quite 
generally within the Ginzburg-Landau theory. Let the elec- 
trodes be described by the (macroscopic) order parameters 
A 1 2 = |A|e' Vl . 2 . Then the phase-dependent part of the energy 
of the junction is E = a[A { Aj + H.c/J + ^IXA^^ + H.c.] 

H where a, 6, . . . depend weakly on T. ]0 Thus for T 

close to T c we estimate 7}Oc|A| 2 oc(r c - T) and 7 2 oc|A| 4 
^{T c — T) 2 , leading to \l 2 /li\<\. With increasing devia- 
tions from ideal geometry \I 2 fI\\ decreases. For large 
enough deviations, negligible values of |/ 2 | are expected. 
These expectations are qualitatively consistent with the ex- 
perimental data (see also Fig. 4). 

So far, our discussion was based solely on symmetry ar- 
guments. Let us attempt a more quantitative analysis of our 
data now. Two different microscopic pictures of asymmetric 
45° Josephson junctions between */-wave superconductors 
have been considered in the literature. The first picture as- 
sumes a microscopically tetragonal material and an ideally 
flat interface. 1 0-1 2 Within this picture, only sample No. 1 can 
be analyzed. Sample No. 2 had I 0 (T= 1.5 K) = 10" 2 M- At 
this temperature only the ^-periodic component of I{<p) was 
observed. At higher temperatures f 0 was not measurable. 
I((p) for sample No. 1 calculated according to the model of 
Ref. 1 1 is shown in Fig. 5. The experimental data can be 
fitted within a relatively broad range of barrier heights. How- 
ever, if we require the I(<p) relation of the 36° junction to be 
fitted by the same (or smaller) barrier height as for the 45° 
junction, we conclude the barrier of the 45° junction to be 
rather low. 14 The dependence of I(<p) on T requires a choice 
of T c ^60 K in the non-self-consistent theory of Ref 11. The 
reduction from the bulk T c = 90 K is probably due to a com- 
bined effect of surface degradation and order-parameter sup- 
pression at the sample surface. The temperature dependence 
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of the ratio of the tt and 277 periodic components in I((p) is 
seen to be in qualitative agreement with experimental data in 
Fig. 3(a). This is explicitly demonstrated in Fig. 4 where we 
compare the experimentally obtained I X 2 with the results of 
the Fourier analysis of the curves in Fig. 5. The divergence 
of 1 2 as T— >0 is an artifact of the ideal junction geometry 
assumed in Ref. 1 1 . If a finite roughness of the interface is 
taken into account, this divergence is cut off and the experi- 
mental data in Fig. 4 do indeed resemble theoretical predic- 
tions for a rough interface. 12 However, the non-self- 
consistent theory of Ref. 11 is unable to explain the 
experimentally observed steep CPR close to the minima of 
the junction energy [see Fig. 3(a)]. 

In a different approach a heavily meandering interface 
with &j= 0j(x) is assumed. Now, the critical current density 
j c (x) is a random function with a typical amplitude 
(\jc(. x )\)~~jc- ^ tne average critical current along the junc- 
tion (j c )<j c , a remarkable 77-periodic component is present 
in the CPR. The relation \I 2 /I\\ depends on the distribution 
of j c (x) and can be much larger than one for (jc)^jc- 15 ' 16 
This model qualitative explains the obtained results for all 
samples, however for a quantitative comparison the actual 
microscopic distribution j c (x) should be known. Note that 



also within the picture of Refs. 15 and 16 the rf-wave sym- 
metry of the pairing state is crucial, otherwise the condition 
Uc)^jc is difficult to satisfy. 

Our present understanding of I(<p) in the asymmetric 45° 
junction is only qualitative. We cannot say whether the re- 
markable amplitude of the 7r-periodic component of /(<p) is 
dominated by the microscopically flat regions, 13 or due to the 
spatial inhomogeneity of the junction. This issue requires 
further study. 

In conclusion, we have measured the magnetic-field re- 
sponse of a single-junction interferometer based on asym- 
metric 45° grain-boundary junctions in YBa 2 Cu 3 0 7 _ x thin 
films, A large 7r-periodic component of 7(<p) has been ex- 
perimentally found, which is in agreement with theoretical 
predictions for c/ JC 2_ v 2-wave superconductors. Hence our re- 
sults provide an additional source of evidence for the t/-wave 
symmetry of the pairing state in the cuprates. 
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For various configurations of Josephson junctions incorporating superconductors with unconventional 
order parameter symmetry, such as most high-TV cuprates, deviations from the standard sinusoidal 
current-phase dependence have been predicted. To this point, these deviations have never been observed 
experimentally We have measured the current-phase relation of high-7 c Josephson junctions, namely, 
YBa 2 Cu 3 0 7 -, thin film bicrystals, comprising symmetric 45° [001] tilt grain boundaries. The current- 
phase relations of all junctions investigated were found to be extremely nonharmomc, in agreement with 
a d* ,2-wave dominated symmetry of the order parameter. [S003 1-9007(98)06674-5] 



PACS numbers: 74.50. +r 

The current-phase relation (CPR) f(<p) is a character- 
istic property of any weak link connecting two supercon- 
ductors. It describes the dependence of the Cooper-pair 
current I p on the phase difference cp of the order parame- 
ters of both superconducting electrodes. In a general form 
it is expressed as 



I P = hf(<p\ 



-1 </(<?) < 1, 



(1) 



I c being the critical current of the weak link. It was 
shown by Josephson [1] that for ideal tunnel junctions 
between conventional superconductors the CPR is sinu- 
soidal, i.e., f(<p) = sin(^). This sinusoidal dependence 
has been confirmed experimentally numerous times for 
standard tunnel junctions between conventional supercon- 
ductors [2]. 

Recently, it has been revealed that the order parameter 
of most high-r c cuprates is unconventional, dominated 
by a d x i- y 2 symmetry component [3-5]. Because of the 
sign change of the order parameter associated with this 
symmetry, strong deviations from the standard sinusoidal 
dependence have been predicted for the current-phase 
relations of various configurations of Josephson junctions 
employing such unconventional superconductors [6-9]. 
In particular, nonharmonic and double-periodic current- 
phase relations are expected for junctions oriented nomi- 
nally perpendicular to the <1 10) direction of one or of 
both electrodes, as well as for junctions for which the 
(110) direction of one of the electrodes is aligned with 
the (100) direction of the other, such as for 45° [001] tilt 
grain boundaries. 

These predictions are highly unusual. Therefore, an ex- 
perimental clarification of the CPR for high-T c junctions 
for which deviations from a standard harmonic behavior 
are expected is desirable. Such experiments will further 



enhance the understanding of the influence of the order 
parameter symmetry on the properties of grain boundaries 
and high-7 c Josephson junctions. In addition, they will 
provide valuable information for the design and use 
of Josephson junction-based circuits, of which many 
characteristics directly depend on the CPR. However, to 
our knowledge, such experiments have not been carried 
out. All available data refer to Josephson junctions 
for which nominally sinusoidal current-phase relations 
are expected. The CPR was measured for weak links 
prepared by ion irradiation [10], for step-edge junctions 
[11,12], and also for 24° bicrystal grain boundaries 
[12]. In nearly all of these cases sinusoidal current- 
phase relations were found. Deviations from a sinusoidal 
dependence have been observed only for one step-edge 
junction, measured at 77 K [12]. These deviations can be 
explained by the influence of thermal noise [13]. 

For these reasons we have investigated the CPR of 
YBa 2 Cu 3 0 7 -^ thin film bicrystals with symmetric 45° 
[001]-tilt grain boundaries, as sketched in Fig. 1(a). For 
these junctions, strong deviations from a sinusoidal CPR 
are anticipated. 

Following a standard approach [14], the CPR was mea- 
sured using a single-junction interferometer configuration 
in which the Josephson junction is part of a superconduct- 
ing loop with a small inductance L. The phase difference 
<p across the junction is controlled by applying an external 
magnetic flux <t> e penetrating the loop: 



<Pe ~ Pfi<p) + <Pn + 27rm . 



(2) 



Here, <p c = 2tt<S? q /^o is the external flux normalized to 
the flux quantum <D 0 (« 2.07 X 10" 15 Tm 2 ). The con- 
stant /3 = 2ttL/ c /<I>o is the normalized critical current, 
(p n is a term accounting for the effective noise, and m is 
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FIG. 1. (a) Schematic representation of a symmetric 45° 
[001]-tilt grain boundary junction in a d x 2- y i superconductor. 
The boundary in the superconducting thin film is meandering, 
leading to the occurrence of tt facets, (b) Schematic of the 
measurement setup. The Josephson junction is denoted by JJ, 
and C indicates the capacitance of the tank circuit. The other 
denotations are explained in the text. 



an integer. As capacitive contributions to the loop current 
are insignificant at the measurement frequencies, the junc- 
tion capacitance has been neglected. The quasiparticle 
current is also negligible, for reasons discussed below. 

The superconducting loop is inductively coupled to a 
tank circuit with inductance Lt [see Fig. 1(b)]. This tank 
circuit is driven with a current I T f at a frequency a> and 
a dc current / dc . Thus, <p e can be expressed as a sum 
of a dc and an rf component <p c — <p<\ c + <p T f. In this 
arrangement, the effective impedance Z e ff (co) of the loop- 
tank circuit combination is a function of cp t . As shown 
by Rifkin and Deaver [14], the CPR can be obtained 
from this dependence, provided that ^? r f <C 1. To obtain 
the CPR for the complete phase range 0 ^ <p ^ 27r, the 
condition /? < 1 has to be fulfilled in addition. 

To enhance the accuracy of the measurement, we have 
adapted this common approach and retrieved the CPR 
from a measurement of the <p& c dependence of the phase 
angle a between the drive current / r f and the tank voltage 
V at the resonant frequency of the tank circuit <dq. As 
described in Ref. [12], at w 0 > the a(<p^ c ) dependence is 
related to the derivative of the CPR /'(<?) ss df((p)/d<p 
in the following way: 



tana(^ dc ) = rr . (3) 

1 + Pf(<p{<Pd c )) 

Here k is the coupling factor between the tank induc- 
tance and the interferometer, k 2 = M 2 /(LLt)> where M 
is the mutual inductance [Fig. 1(b)]. Using Eq. (3), from 
the measured a(<Pdc) dependence f f (<p(<Pdc)) is obtained. 
The CPR is restored by integrating f f (<p(<Pd C )) numeri- 
cally, using the d(p{(p^ c )/d<p^ c dependence obtained from 
differentiating Eq. (2) with respect to <p^ c . 

The samples investigated consisted of three bicrys- 
talline YBa 2 Cu 3 0 7 -* films with a T c (R = 0) of 88 K. 
The films, with thickness t = 100 nm, were deposited by 
standard pulsed laser deposition on (OOl)-oriented SrTiC>3 
bicrystalline substrates [15] with symmetric [001]-tilt mis- 
orientation angles of 45° ± 2° and were subsequently 
patterned by Ar ion-beam etching into 8X8 mm 2 or 
5X5 mm 2 square washer single-junction interferome- 
ter structures. The widths of the junctions were b 
2-3 jim. The washer holes had a side length of 50 /xm, 
leading to L « 80 pH. 

To minimize the influence of external noise, the 
samples were measured in superconducting and double 
magnetic shielding at a temperature of 4.2 K. The 
condition {$ < 1 for the investigated interferometers 
was confirmed experimentally from the character of its 
response versus (p^ c [12], 

YBa2Cu307-^ grain boundaries with a symmetric [001] 
tilt angle of 45° typically have a normal-state interface- 
resistivity p n > 1 X 10~ 8 ft cm 2 , which we also mea- 
sured for boundaries fabricated under identical conditions 
as the junctions used in the present experiments. In the 
configuration used, this p„ corresponds to normal-state 
resistances R n > 1 ft. Accordingly, the relaxation time 
of the interferometer r = L/R is short (r l/o>o), and 
hence the quasiparticle current is negligible [12]. 

For the measurements of a(<p^ c ), two tank circuits with 
quality factor Q = 120 and inductance Lt = 0.4 /iH, 
qjq = 30 MHz, and L T = 0.73 jjlU, oj 0 = 23 MHz, re- 
spectively, were employed. The phase angle was recorded 
as a function of /dc after amplification of the tank voltage 
by a high-impedance amplifier. The coupling coefficient 
k was determined from the period of the a(/d C ) depen- 
dence. Values of k = 0.072 and k = 0.054 for the re- 
spective tank circuits are obtained. To ensure the validity 
of the small-signal limit, the measurements were carried 
out with <p x f < 0.15. 

A typical a(/dc) dependence is shown in Fig. 2. The 
corresponding CPR, depicted in Fig. 3, is clearly deviat- 
ing from the standard sinusoidal behavior. Samples fabri- 
cated on different substrates and measured with both tank 
coils followed closely the same behavior. It is empha- 
sized that the experimental setup employed and the pro- 
cedure followed are identical to those used to measure the 
current-phase relations of step-edge junctions and thin- 
film bicrystals with a symmetric [001] tilt of 24°. For 
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FIG. 2. Phase angle a between the driving current and the 
output voltage measured at 4.2 K as a function of the dc current 
/do for an YBa2Cu307_* single junction interferometer circuit 
containing a symmetric 45° [001]-tilt grain boundary. 



all of those samples nominally sinusoidal current-phase 
relations were observed, and all deviations of the apparent 
CPR from a sinusoidal one can be attributed to thermal 
noise [12,13]. 

The measured deviations from a sinusoidal de- 
pendence for the current-phase relations of these 
45 ° bicrystals are startling. It is important to 
note that the effective Josephson penetration depth 
Aj = [ < I ) o/(47r/to(7c)A)] 1 / 2 » 5 jxm is larger than the 
width of the junction b (narrow-junction limit). Here 
A is the London penetration depth. Although several 
mechanisms are known to cause nonsinusoidal current- 
phase relations for narrow junctions fabricated from 
conventional superconductors, all of these mechanisms 
fail to account for the anomalous dependencies presented. 




-l.a 



FIG. 3. The normalized current through the junction (3f{<p) 
as a function of the phase difference <p restored from the 
measured a(/d C ) as shown in Fig. 2. For comparison, the 
function f$ sin(^>) with /3 = 0.8 is plotted as a solid line. 



First, one potential source of such deviations is thermal 
noise. To evaluate its influence we consider a sinusoidal 
CPR and calculate with Eq. (3) the a(<p^ c ) dependence, 
assuming a thermally induced Gaussian spread p(<p n )- 
With this, the value of tana(<pdc) is given by [12] 

Xma(<p dc ) = k 2 Q/3 



/CO 



cos <p(<p t ,<p n ) 

1 + f3 COS <p(<Pe,<p n ) 



p{<Pn)d<p n , 



(4) 

Using Eq. (4), minor deviations of the CPR from a stan- 
dard sinusoidal behavior are well described quantitatively 
for 24° boundaries measured at 77 K [12]. However, no 
realistic set of f3 and <p n exists to account for the large 
deviations of the CPR observed for the 45° boundaries. 

Second, for weak links with a high current density, 
current-induced suppression of the order parameter in 
the electrodes close to the weak link can give rise to 
nonharmonic current-phase relations [16]. It is unrealistic 
that this effect is the cause for the deviations presented 
here, as the intragrain critical current density exceeds 
10 7 A/cm 2 at 4.2 K and is therefore much larger than 
the grain boundary (j c ) < 4000 A/cm 2 at the same 
temperature. 

Third, several additional mechanisms, described in [2], 
lead to deviations from a harmonic CPR. In all these 
cases, the slope of f(<p) at <p = 0 is smaller than at <p = 
77, which is in contrast to our results. Therefore, these 
mechanisms cannot explain the current-phase relations 
observed either. 

On the other hand, as will be pointed out in the 
following, the measured CPR can be accounted for 
by the unconventional order parameter symmetry of 
YBa2Cu307-jc and by the microstructural properties of 
the grain boundaries, in particular by their faceted nature 
[17,18]. Interestingly, due to the d x 2- y 2 -wave character 
of the order parameter, the faceting has a more significant 
influence on the electronic properties of boundaries with 
a misorientation close to 45° than on boundaries with 
considerably smaller misorientation angles [18]. This 
also concerns the CPR, as the 45° boundaries contain a 
higher density of facets that themselves show anomalous 
behavior. Two kinds of such anomalies, to be considered 
here, are described in the literature. 

First, due to the sign difference of the adjacent lobes of 
the d^-^ -wave order parameter, many facets are biased 
with an additional it phase shift (it facets) [17-19]. 
These phase shifts give rise to unconventional junction 
properties, such as a spontaneous generation of magnetic 
flux in the grain boundary junction [19-21]. As described 
in Ref. [21], the local phase difference <p(x) along the 
grain boundary (0 < x < b) can be written as 

+ (5) 

where f (jc) is a rapidly alternating function accounting 
for the 7r phase shifts and for the spontaneously generated 
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magnetic flux in the junction, and ft(x) is the remaining 
slowly varying phase difference reflecting the magnetic 
flux in the interferometer loop. For a narrow junction 
(b < A;), ft is independent of x. With Eq. (5), the time- 
independent sine-Gordon equation, describing the spatial 
dependence of the local phase difference over the junc- 
tion, becomes 

2 d 2 £{x) j c (x) . r , w Yl i£\ 

A; = jjj sm[ft + (6) 

The solution of this equation, and thus the pattern 
of self-generated flux, depends on ft. The redistribution 
of this flux by a change of ft is expected to lead to 
remarkable deviations from a harmonic dependence for 
the CPR measured for the entire junction, also if the local 
CPR is nominally sinusoidal. 

Second, it has been proposed [7,8] that the CPR of 
facets formed by the (110) and by the (100) planes 
of the adjacent grains is periodic with tt and thus has 
a double periodicity as compared to the standard case. 
Transmission electron microscopy has revealed that 45° 
[001] tilt grain boundaries in YBa2Cu307-jc tend to be 
composed for a considerable part of such facets [22]. For 
the whole junction, this leads to an anomalous CPR: 

/ = I c \ sin ft + I C 2 sinlft , (7) 

by which the observed CPR can be described. 

In summary, the current-phase relations of grain bound- 
ary junctions with a misorientation of 45° were mea- 
sured with a modified Rifkin-Deaver method. The CPRs 
were deduced from measurements of the phase angle be- 
tween the rf drive current and the rf tank voltage. The 
current-phase relations of the Josephson junctions showed 
pronounced deviations from a harmonic behavior, which 
cannot be accounted for by thermal noise or by other stan- 
dard mechanisms, but are attributed to the d^-^-wave 
symmetry of the order parameter and the faceting of the 
grain boundaries. 
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The predominant d-wave pairing symmetry in high-temperature superconductors allows for a 
variety of current-phase relations in Josephson junctions, which is to a certain degree fabrication 
controlled In this Letter, we report on direct experimental observations of the effects of a non- 
sinusoidal current-phase dependence in YBCO dc SQUIDs, which agree with the theoretical description 
of the system. 



DOI: 10. 1103/PhysRevLetL90. 117002 

It is well established [1] that the wave function of a 
Cooper pair in most cuprate high-temperature supercon- 
ductors (HTS) has a </-wave symmetry. Its qualitative 
distinction from, e.g., the anisotropic s-wave case is that 
the order parameter changes sign in certain directions, 
which can be interpreted as an intrinsic difference in the 
superconducting phase between the lobes equal to 7r. 

The latter leads to a plethora of effects, such as forma- 
tion of Andreev bound states at surfaces and interfaces in 
certain crystallographic orientations [2-4]. The current- 
phase dependence I s (<f>) in Josephson junctions formed by 
dd junctions, as well as by sd junctions comprised of a 
cuprate and a conventional superconductor, depends both 
on the spatial orientation of the d-wave order parameter 
with respect to the interface, and on the quality of the 
latter [5-9]. Time-reversal symmetry can also be sponta- 
neously violated and thus spontaneous currents generated 
[10-12]. Another effect can be doubling of the Josephson 
frequency [6,13,14]. 

In this Letter, we report on experimental observa- 
tions of strong effects of an unconventional current-phase 
relation on the dynamics of two dd junctions integrated 
into a superconducting interference device (SQUID) 
configuration. 

Since /$(<£) must be a 27r-periodic odd function, it can 
be expanded in a Fourier series. In most cases, only the 
first two harmonics give a significant contribution to the 
current: 

h(4>) = 7 c sin^ - /? sin20. (1) 

In Josephson systems of conventional superconductors, 
the second harmonic will usually be negligible [15] but 
in dd junctions the second harmonic may dominate. If 
/J 1 > /J/2, the equilibrium state is no longer <f> = 0 but 
becomes double degenerate at <f> = ± arccos(/J/2/J?) — 
±77-/2. The system can then spontaneously break time- 
reversal symmetry by choosing either state. Spon- 
taneous currents as well as fluxes can be generated in 
this state. The potential will have the shape of a double 



PACS numbers: 74.50.+r, 85.25.Dq 

well, and there are reasons to believe that it will be 
possible to observe quantum coherence in this system. 
The presence of a second harmonic in the current-phase 
relation (CPR) of a dd junction was confirmed by 
Il'ichev et al [8]. 

A nonsinusoidal CPR of the junctions will change the 
dynamics of a dc SQUID [16]. Regarding the junctions as 
magnetically small, the supercurrent through the SQUID 
in the presence of an external flux & x = 3> 0 * (^jc/2w) 
can be written as 

hit* <t> x )-i\A sin^-/? lS in(2*) 

+ I\ 2 sin(<£ + <f> x ) - /° 2 sin2(0 + </> x ). (2) 

The critical current through the SQUID is given by the 
usual expression l c {<i> x ) = max^/,^, <f> x ). The time- 
averaged voltage over the SQUID in the resistive regime 
is readily obtained in the resistively shunted junction 
approximation. By introducing 8[<f>, <f> x ] = <f> 2 ~ <f>\ and 
applying the same method as in [17] with the necessary 
generalizations, we obtain the following for the average 
voltage over the SQUID: 

-'■K)-"H)r 

(3) 

Here G\ 2 are the normal conductances of the junctions, 
and 

8 + <t> x + ^ [/ 2 (<£ - 5/2) - /, {<f> + 8/2)] = 0 (4) 

gives the difference, 3, in phase drops across each junc- 
tion. In deriving (3) and (4), we have assumed that the 
inductance L is equally divided between the SQUID 
arms. We have also neglected the spontaneous magnetic 
fluxes in the dd junctions, due to their small amplitude 
[11,18]. Though (4) is only explicitly solvable in the limit 



117002-1 0031-9007/03/90(11)/117002(4)$20.00 



© 2003 The American Physical Society 



117002-1 



PHYSICAL REVIEW LETTERS 
Volume 90, Number 11 . — 



L — ► 0, it always yields -<f> x ] = -S[<f> t <f> x l This 

means that the usual inversion symmetry is retained. 

The results of numerical calculations based on (2) and 
(3) are shown in Fig. 1. The cusps in the critical current 
correspond to the points at which the global maximum in 
(1) switches from one local maximum to another [16]. 
Note the quasi-<& 0 / 2 periodicity of the current isolines in 
the V - <l> x picture, reflecting the current-phase depen- 
dence (1), and their shift along the <& x axis, which de- 
pends on the sign of the bias current (as it must to 
maintain the central symmetry with respect to the ori- 
gin). The shift does not depend on the magnitude of the 
current since we neglect the self-inductance. For large 
biases, the <& 0 periodicity is restored. Indeed, as the bias 
grows, one set of minima of the washboard potential, 
U = (h/2e)[-I l co&4> + (/ I! /2)cos2<£ - /<£], disappears 
first unless the first harmonic I 1 is exactly zero. 

We have fabricated and studied a large number of dc 
SQUIDs. The samples were fabricated from 250 nm thick 
YBCO films deposited on SrTi0 3 bicrystals. The grain- 
boundary junctions (GBJs) are of the asymmetric 
[001]-tilt type with the misorientation angle of 45° 
(0°-45° GBJ). For more information on GBJs see, for 
example, Ref. [19]. 

The pattern was defined using E-beam lithography and 
then transferred to a carbon mask employing a multistep 
process. Finally, the YBCO is etched through the mask 
using ion milling. This scheme allows us to fabricate 




'±2 -1 5 -1 -05 0 0.5 1 1 5 2 

T x 0 

FIG. 1. The results of simulations of the I c - </> x and V - <f> x 
dependence for a dc SQUID with l\ x = 1, 4 = 0.1, = 0.2 
and l\ = 0.4 (arbitrary units). The different curves correspond 
to bias currents in the range / = l\ x to / - 5l\ v We assume 
L — 0 and G\ = G 2 - 
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high-quality bicrystal junctions as narrow as 0.2 /xm, as 
has been reported elsewhere [20]. In the SQUIDs under 
investigation, the junctions are nominally 2 fim wide; 
hence, the fabrication-induced damage of the junctions 
is small. 

The measurements were done in an EMC-protected 
environment using a magnetically shielded LHe cryostat 
However, the magnetic shielding is imperfect, as is evi- 
dent from the fact that the expected zero-field response of 
our SQUIDs is not exactly at zero. The measuring elec- 
tronics is carefully filtered and battery powered whenever 
possible. In order to measure the dependence of the 
critical current on the applied field, we used a voltage 
discriminator combined with a sample-and-hold circuit 
All measurements reported here were performed at 4.2 K. 

The SQUID loops are (15 X 15) /un 2 . The numerically 
calculated inductance [21] is approximately 25 pH, yield- 
ing the factor (3 = IttLIJ^o between 0.5-2. 

The SQUIDs were largely nonhysteretic with a resis- 
tance of about 2 fl. The measured critical current varies 
from sample to sample but is in the range of tens of 
microamperes giving a current density of the order of 
j = iq3 A/cm 2 . The estimated Josephson penetration 
length A 7 - 4> 0 />/Wo)cAl is approximately 2 fim in 
all junctions, which means that the junctions are magneti- 
cally short This is supported by the quasiperiod of the 
pattern in Fig. 2 being close to the expected value 
<f> 0 /2A L w [17]. The differential conductance curves do 
not show any trace of a zero bias anomaly (ZBA), as is 
expected for 0°-45° GBJs. ZBAs have been observed by 
other groups in GBJs with other orientations [2]. 

The critical current is plotted as a function of applied 
magnetic field for two SQUIDs in Fig. 3. The result is in 
qualitative agreement with theory if we assume that the 
SQUID junctions have different ratios of the first and 
second harmonics of the critical current This assumption 




Applied Magnetic Field (ml) 



FIG. 2. Critical current as a function of magnetic field at 
4.2 K. The dashed box indicates the area plotted in Fig. 3(a). 
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is supported by the fairly small modulation depth [it is 
easy to see from Eq; (2) that I c would go exactly to zero in 
a SQUID with junctions of identical / c2 // c |]* 

We can fit the data to Eq. (2), if we compensate for the 
residual background magnetic field and assume that we 
have a small excess current (of the order of a few /xA) in 
the junctions. The fitting parameters again confirm that 
there is a large asymmetry between the arms of the 
SQUIDs. Note that the model does not consider the flux 
penetration into the junctions, 

The result for fields of the order of mT is presented in 
Fig. 2, which shows the I c modulation of the SQUID 
enveloped by an anomalous Fraunhofer pattern quite 
similar to what has been reported by other groups 
[22,23] for 0°^5° GBJs. Note the inversion symmetry 
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FIG. 3. Critical current as a function of applied magnetic 
field for two different SQUIDs that are nominally identical. 
The solid line represents the fitted expression. The fitting 
parameters are as follows: (a) l[ x — 9 /xA, l\ 2 = 0.3 /*A, l\ = 
3.7 M, and 7° = 22.7 /*A; (b) - 7.8 /xA, 7* 2 - 3.0 /xA, 
/iij =53 m a, and /° 2 = 4.3 ^A. In both cases, the fit has been 
adjusted with respect to the residual background field and the 
excess current of the junctions. 



of the pattern with respect to the origin. That the global 
maximum is not in the center can be explained in several 
ways; it has been shown, for example, that this could be 
due to the presence of so-called 7r loops in the junction 
interface [24]. 

Figure 4 shows the V-B dependence of one of the 
SQUIDs. The pattern is again field inversion symmetric. 
The overall structure is the same as in the model depen- 
dence of Fig. 1, but there is also an additional shift due to 
self-field effects, which depends on the magnitude of the 
bias current and corresponds (at maximum) to a flux 
~0.14> 0 . In a beautiful experiment, a similar dependence 
was recently observed by Baselmans et al in a Nb-Ag- 
Nb SNS junction where current injectors were used to 
change the occupation of current-carrying states in the 
normal region [25]. A deviation from the model occurs at 
V = 100 /tV where the minima and maxima switch. This 
is probably due to an LC resonance in the SQUID. Taking 
L = 25 pK, this would require C = 0.8 pF, which agrees 
with our measurements on single junctions 

Remarkably, the observed offset of the V-B character- 
istics with respect to the direction of the bias current 
appears to be a much more robust manifestation of the 
presence of a second harmonic of the Josephson current 
than the shape of the I c - B curves itself. We observed the 
shift even in SQUIDs with the smallest junctions down to 
0.5 /xm wide, where the deviations from the usual sinu- 
soidal CPR were not obvious from the I c - B dependence. 

Generally, the nature of the transport through a GBJ 
will depend on its transmissivity D. Il'ichev et al [8] 
have reported values of D as high as 0.3 in symmet- 
ric (22.5°-22.5°) dd junctions as opposed to the usual 
estimate for a GBJ, D - 10" 5 -10~ 2 . Since usually 
/? jl\ a D, a high-transmissivity GBJ is required in order 




-10 -5 0 5 

Applied Magnetic Field (^iT) 

FIG. 4. Voltage modulation as a function of applied magnetic 
field for the SQUID whose / c - B is shown in Fig. 3(a). The 
pattern is again inversion symmetric. Note the sign change at 
100 /tV, which we believe is due to a LC resonance in the 
SQUID loop. 
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to observe effects of the second harmonic. An estimate of 
the average transmissivity of our junctions would be 
p ab l/R N A ~ 10~ 2 [26] assuming /, the mean-free path, 
to* be equal to 10 nm and a resistivity in the a-b plane p ab 
equal to 10" 4 Hem. This is still too low to explain the 
strong second harmonic we observe. However, it is known 
from, e.g., TEM studies [19], that the grain-boundary is 
far from uniform; the properties can significantly vary 
depending on the local properties of the interface, effects 
such as oxygen diffusion out of the GB, etc., which are 
difficult to control. It is therefore reasonable to assume 
that there are many parallel transport channels through 
the GB [27,28]. Channels with high-transmissivity domi- 
nate the transport and might have D ~ 0.1 even though 
the average transmissivity is much lower. This is also 
consistent with the fact that most of our SQUIDs seem 
to be highly asymmetrical which is to be expected if the 
distribution of channels is random. The ratios of l\ and /? 
can vary as much as 10 times between two junctions in the 
same SQUID, even though the fluctuations of the total I c 
from sample to sample are much smaller. It is also clear 
from general considerations that a high value of /? ex- 
cludes a high value of l\, since the second harmonic 
usually dominates if the odd harmonics of the super- 
current are canceled by symmetry [29]. 

Recent studies of 0°-45° GBJs have demonstrated that 
the SQUID dynamics can be altered by the rf-wave order 
parameter in YBCO [30]. It is, however, important to 
point out that our results do not directly relate to, e.g., 
tetracrystal tt-SQUID experiments; the latter crucially 
depend on having one ir junction with negative critical 
current, but still only the first harmonic present in I c {<f>). 
Our SQUIDs have a conventional geometry, but uncon- 
ventional current-phase relations. 

One explanation for the pronounced effects of the 
second harmonic could be that relatively large sections 
of the interface are highly transparent and have a low de- 
gree of disorder. This in turn could be related to our fabri- 
cation scheme which seems to preserve the integrity of 
the barrier. This makes feasible their applicability in the 
quantum regime and supports our expectations that quan- 
tum coherence can be observed in this kind of structures. 

In summary, we have observed a very pronounced sec- 
ond harmonic in the current-phase relation of a "conven- 
tional" YBCO dc SQUID with 0°-45° grain-boundary 
junctions. It has strongly influenced the SQUID dynam- 
ics. All details of the SQUID behavior were explained 
within a simple model of a dd junction with relatively 
high transparency. We believe that these effects are im- 
portant for better understanding of HTS Josephson junc- 
tion and SQUIDs. 
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We present various concepts and experimental procedures to produce biepitaxial YBa2Cu3C>7-x 
grain boundary Josephson junctions. The device properties have an interesting phenomenology, 
related in part to the possible influence of tl 7r- loops". The performance of our junctions and Su- 
perconducting Quantum Interference Devices indicates significant improvement in the biepitaxial 
technique. Further, we propose methods for fabricating circuits in which "0-" and "7r-loops" are 
controllably located on the same chip. 



I. INTRODUCTION 

The possibility of realizing electronic circuits in which 
the phase differences of selected Josephson junctions are 
biased by w in equilibrium is quite stimulating. 1 The 
concept of such 7r-phase shifts was originally developed 
in the "extrinsic" case for junctions with ferromagnetic 
barriers 2 and in the "intrinsic" case for junctions exploit- 
ing superconductors with unconventional order parame- 
ter symmetries. 3 As a result of the possible d x 2_ y 2 or- 
der parameter symmetry of high critical temperature su- 
perconductors (HTS), 4 the presence of intrinsic it loops 
has also been considered for HTS systems. 5 This has 
been discussed recently in view of novel device con- 
cepts, and in particular for the implementation of a 
solid state qubit 1,6-8 and for Complementary Joseph- 
son junction electronics. 9 In this paper we discuss how 
YBa2Cu307_a: (YBCO) structures made by the biepi- 
taxial technique 10,11 can be successfully employed to pro- 
duce arbitrary circuit geometries in which both "0" and 
7r-loops are present, and possibly to obtain a doubly de- 
generate state. 1,6 Of course, great caution should be used 
because of stringent requirements on junctions parame- 
ters for practical applications of such devices. 

Josephson junctions based on artificially controlled 
grain boundaries have been widely employed for funda- 
mental studies on the nature of HTS. 4 ' 7,8 The lack of a 



reliable technology based on the traditional trilayer con- 
figuration (i.e. a sandwich type junction with an insu- 
lator between the two superconducting electrodes) also 
enhanced interest in GB Josephson junctions for appli- 
cations. Although the mechanism of high-Tc supercon- 
ductivity and the influence of grain boundaries on the 
transport properties are not completely determined, re- 
producible and good quality devices are routinely fab- 
ricated. YBCO GB junctions are usually classified as 
bicrystals, 12 biepitaxials, 11 and step-edges, 13 depending 
on the fabrication procedure. The bicrystal technique 
typically offers junctions with better performances and 
allows in principle the realization of all different types 
of GBs ranging from [001] and [100] tilt to [100] twist 
boundaries. GB junctions based on the step-edge and 
biepitaxial techniques offer the advantage, with respect 
to the bicrystal technology, of placing the junctions on 
the substrate without imposing any restrictions on the 
geometry. A comparison between the different GB tech- 
niques is far beyond the aim of this paper. Nevertheless 
we intend to show that significant improvements with re- 
spect to the original technique developed by Char et al. 11 
are possible for biepitaxial junctions, and that the result- 
ing devices have potential for applications. As a matter 
of fact, in traditional biepitaxial junctions, the seed layer 
used to modify the YBCO crystal orientation on part of 
the substrate produces an artificial 45° [001] tilt (c-axis 



tilt) GB. The nature of such a GB seems to be an intrinsic 
limit for some real applications. A convincing explana- 
tion has been given in terms of the d-wave nature of the 
order parameter and more specifically by the presence 
of 7r-loops. 14 As demonstrated by studies on bicrystals, 
based on the same type of 45° [001] tilt GB, the presence 
of 7r -loops reduces the Ic^n values (where Ic and Rjv 
are the critical current and the high normal state resis- 
tance respectively), produces a dependence of the critical 
current Ic on the magnetic field H quite different from 
the Praunhofer-like pattern, and generates unquantized 
flux noise at the grain boundary. 14 

We will show that the implementation of the biepi- 
taxial technique 10 we developed to obtain 45° [100] tilt 
and twist (a-axis tilt and twist) GBs junctions makes 
such a technique interesting for both applications and 
fundamental studies. The phenomenology observed for 
the junctions based on these GBs and Scanning SQUID 
Microscopy investigations demonstrate the absence of ir- 
loops, as we expect from their microstructure. As a con- 
sequence higher values of the I c Rn values, a FVaunhofer 
like dependence of Ic on the magnetic field and lower 
values of the low frequency flux noise, when compared 
with 45° c-axis tilt GBs, have been measured. These 
features are important tests to employ junctions for ap- 
plications. Scanning SQUID Microscopy investigations 
also gave evidence of "fractional" vortices in the pres- 
ence of impurities. Finally, we extended the biepitaxial 
process to other types of GB by using different seed layers 
to obtain junction configurations where 7r loops can be 
controllably produced. We shall not dwell on conceptual 
principles and actual feasibility of qubit devices. Instead 
we discuss the importance of the biepitaxial technique in 
having "0" and "tt" loops on the same chip. This makes 
the biepitaxial technique more versatile and promising 
for circuit design. 



of the GB interface, we have also successfully employed 
vicinal substrates. However, most of the transport prop- 
erties presented in this paper refer to samples not using 
vicinal substrates. Detailed structural investigations on 
these GBs, including Transmission Electron Microscopy 
(TEM) analyses, have been performed and the results 
have been presented elsewhere. 10,16 

Depending on the patterning of the seed layer and the 
YBCO thin film, different types of GBs ranging from the 
two ideal limiting cases of 45° a-axis tilt and 45° a-axis 
twist have been obtained (see Fig.l). The intermediate 
situation occurs when the junction interface is tilted at 
an angle a different from 0 or 7r/2 with respect to the a- 
or b-axis of the [001] YBCO thin film. In all cases, the or- 
der parameter orientations do not produce an additional 
7T phase shift along our junction, in contrast with the 
45° asymmetric [001] tilt junctions. As a consequence, 
no 7r loops snould occur independently of the details of 
the interface orientation. In Fig. 1 we consider ideal in- 
terfaces and neglect meandering of the GBs or interface 
anomalies that will be considered below. The Ce0 2 seed 
layer may produce a more complicated GB structure, in 
which a 45° c-axis tilt accompanies the 45° a-axis tilt or 
twist (see Fig.2a). 15 In this case, as shown in Fig.2b, 7r 
loops should occur in analogy with the traditional biepi- 
taxial junctions based on 45° c-axis tilt GBs. In both 
Figs. 1 and 2 we display the possible d x 2_ y 2 -wave or- 
der parameter symmetry in the junction configuration. 
Junctions were typically 4 microns wide. We also per- 
formed systematic measurements on SQUIDs based on 
the structure employing MgO as a seed layer and SrTi0 3 
as a substrate. DC SQUIDs in different configurations 
and with loop inductance typically ranging from 10 to 
100 pH have been investigated. The typical loop size 
leading to the 10(100) pH inductance is approximately 
10 2 fim 2 (lOVm 2 ). 



II. DEVICES: CONCEPTS AND FABRICATION 
PROCEDURE 

As mentioned above, the biepitaxial technique allows 
the fabrication of various GBs by growing different seed 
layers and using substrates with different orientations. 
We have used MgO, Ce0 2 and SrTi0 3 as seed layers. The 
MgO and Ce0 2 layers are deposited on (110) SrTi0 3 sub- 
strates, while SrTi0 3 layers are deposited on (110) MgO 
substrates; in all these cases the seed layers grow along 
the [110] direction. Ion milling is used to define the re- 
quired geometry of the seed layer and of the YBCO thin 
film respectively, by means of photoresist masks. YBCO 
films, typically 120nm in thickness, are deposited by in- 
verted cylindrical magnetron sputtering at a temperature 
of 780° C. YBCO grows along the [001] direction on MgO 
(substrates or seed layers) and on the Ce0 2 (seed layers), 
while it grows along the [103]/[013] direction on SrTi0 3 
(substrates or seed layers). In order to select the [103] or 
[T03] growth and to ensure a better structural uniformity 



III. EXPERIMENTAL RESULTS 

A. Biepitaxial junctions employing MgO seed layers 

In this section we attempt to cover most of the phe- 
nomenology of the transport properties of 45° a-axis tilt 
and twist biepitaxial junctions. In Fig. 3, current vs volt- 
age (I-V) characteristics of a typical biepitaxial junction 
are given for various temperatures close to the critical 
temperature. In the inset the corresponding I-V charac- 
teristic at T = 4.2 K is reported. They are closely de- 
scribed by the resist ively-shunted-j unction (RSJ) model 
and no excess current is observed. Nominal critical cur- 
rent densities J c of 5 x 10 2 A/cm 2 at T = 77 K, and of 9 x 
10 3 A/cm 2 at T = 4.2 K have been measured respectively. 
The Rjv value (3.2 Q) is roughly independent of the tem- 
perature for T < Tc, providing a normal state specific 
conductance a N = 70 (/xOcm 2 ) -1 . The maximum work- 
ing temperature Tc of this device was 82 K. In this case 
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l c R N is 1.3 mV at T = 4.2 K. These values typically 
ranged from 1 mV to 2 mV at T = 4.2 K. They are larger 
for the corresponding Jc values than those provided by 
conventional biepitaxials, and are of the same order of 
magnitude as in GB bicrystal and step edge junctions. 10 
While the values of critical current density and normal 
state specific conductance in the tilt case are quite dif- 
ferent from the twist case, the I c Rjv values are approxi- 
mately the same for both. Moreover IcRjv does not scale 
with the critical current density. 10 In the tilt cases Jc ~ 
0.5-10 10 3 A/cm 2 and <r N ~ 1-10 (^cm 2 )" 1 are mea- 
sured at T = 4.2 K respectively. Twist GBs junctions 
are typically characterized by higher values of Jc in the 
range 0.1-4.0 x 10 5 A/cm 2 and of a N in the range 20-120 
(/xHcm 2 )" 1 (at T = 4.2 K). For the twist case deviations 
from the RS J model are more marked as a result of higher 
critical current densities. For high values of Jc GB junc- 
tions do not present any clear modulation of the critical 
current as a function of the magnetic field. 

A demonstration of the possibility of tailoring the criti- 
cal current density and of the different transport regimes 
occurring in the tilt and twist cases has been given by 
measuring the properties of junctions with different ori- 
entations of the GB barrier on the same chip. By pattern- 
ing the seed layer as shown in Fig.4a, we could measure 
the properties of a tilt junction and of junctions whose in- 
terface is tilted in plane by an angle a = 30°, 45° and 60° 
with respect to the a- or b-axis of the [001] YBCO thin 
film respectively. In all cases the order parameter orien- 
tations do not produce an additional it phase shift along 
our junction, in contrast with the 45° [001] tilt junctions, 
and no tt loops should occur. We measured the expected 
increase of the critical current density with increasing an- 
gle, which corresponds on average to an increase of the 
twist current component. The values measured at T = 
4.2 K are reported in Fig. 4a and range from the min- 
imum value Jc = 3 x 10 2 A/cm 2 in the tilt case to the 
maximum Jc = 10 4 A/cm 2 corresponding to an angle of 
60°, for which the twist component is higher. The con- 
sistency of this result has been confirmed by the values of 
normal state resistances, which are higher in the tilt case 
and decrease with increasing a. The IcR/v values are 
about the same for all the junctions independently of the 
angle a. In Fig. 4b the I-V characteristics measured at 
T = 4.2 K, corresponding to the junctions of Fig.4a, are 
shown for approximately the same voltage range. Devia- 
tions from RSJ behavior appear for higher values of the 
critical current density (a = 60°). These results demon- 
strate that the grain boundary acts as a tunable bar- 
rier. This possibility of modifying the GB macroscopic 
interface plane by controlling the orientation of the seed 
layer's edge is somehow equivalent to the degree of free- 
dom offered by bicrystal technology to create symmetric 
or asymmetric GBs, with the advantage of placing all the 
junctions on the same substrate. The 45° a-axis tilt and 
twist GBs and the intermediate situations can represent 
ideal structures to investigate the junction physics in a 
wide range of configurations. The anisotropy of the (103) 



films and the possibility to select the orientation of the 
junction interface by suitably patterning the seed layer, 
and eventually the use of other seed layers which produce 
different YBCO in plane orientations, allow the fabrica- 
tion of different types of junctions and the investigation 
of different aspects of HTS junction phenomenology. In 
particular we refer to the possibility of changing the tun- 
neling matrix elements (by selecting the angle a) and 
to use the anisotropy of the layered structure of YBCO 
properties and of the order parameter symmetry. 

The study of the junction properties in the presence 
of an external magnetic field H is a fundamental tool 
for the investigation of the Josephson effect in the vari- 
ous junctions, as well as a test of junction quality. 17 We 
observe modulations of the critical current Ic following 
the usual Fraunhofer-like dependence. The Ic (H) pat- 
terns are mostly symmetric around zero magnetic field, 
and in all samples the absolute maximum of Ic occurs 
at H=0. The presence of the current maximum at zero 
magnetic field is consistent with the fact that in our junc- 
tion configuration the order parameter orientations do 
not produce an additional 7r phase shift, in contrast with 
the 45° [001] tilt GB junctions. 14 ' 10 Some examples are 
given in Fig. 5, where the magnetic pattern relative to 
a SQUID and a single junction at T = 4.2 K are shown 
respectively. In the former case we can also distinguish a 
smaller field modulation (with a period of 8 mG) which 
corresponds to the SQUID modulation (inset a). In the 
latter case the I-V characteristics are reported for differ- 
ent magnetic fields (inset b). Despite the Fraunhofer-like 
dependence, some deviations are evident, in agreement 
with most of the data available in literature. 

For sake of completeness we also acknowledge some 
work we carried out by investigating Fiske steps as a 
function of H in other junctions, giving some evidence of 
a dielectric-like behavior 18 of some of the layers at the 
junction interface. We already reported about this work 
elsewhere. 19 The Fiske steps do not depend on the use 
of a particular substrate, since they have been observed 
in junctions based both on SrTi0 3 and MgO substrates. 
Typical values of the ratio between the barrier thickness 
t and the relative dielectric constant e r range from 0.2 
nm to 0.7 nm. Considerations on the dependence of Ic 
on the temperature (T) can be also found in Ref. 19 . In 
junctions characterized by lower critical current densities, 
Ic tends to saturate at low temperatures, in contrast to 
those characterized by higher critical currents, for which 
there is a linear increase. 10,20 



B. Scanning SQUID microscopy on biepitaxial 
junctions with MgO seed layer 

Figure 6 is a scanning SQUID microscope 21 image of 
a 200x200jum 2 area along a grain boundary separating a 
(100) region from a (103) region (as labelled in the fig- 
ure) of a thin YBCO biepitaxial film grown as described 
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above. The position of the grain boundary is indicated 
by the dashed line. The image was taken at 4.2K in liquid 
helium with an octagonal SQUID pickup loop 4 microns 
in diameter after cooling the sample in a few tenths of a 
fiT externally applied magnetic field normal to the plane 
of the sample. The grey-scaling in the image corresponds 
to a total variation of 0.13$ 0 of flux through the SQUID 
pickup loop. Visible in this image are elongated inter- 
layer Josephson vortices in the (103) area to the right, 
and "fractional" vortices in the (100) area to the left, of 
the grain boundary. Fits to the interlayer vortices give 
a value for the c-axis penetration depth of about 4jxm. 
The "fractional" vortices are spontaneously generated in 
the (100) film, regardless of the value of external field 
applied. 22 Temperature dependent scanning SQUID mi- 
croscope imaging shows that this spontaneous magneti- 
zation, which appears to be associated with defects in the 
film, arises when the film becomes superconducting. 4,23 
Although it is difficult to assign precise values of total 
flux to the "fractional" vortices, since they are not well 
separated from each other, fits imply that they have less 
than $o of total flux in them, an indication of broken 
time-reversal symmetry. Although there is apparently 
some flux generated in the grain boundary region, the 
fact that these SQUIDs have relatively low noise seems 
to indicate that this flux is well pinned at the tempera- 
tures at which the noise measurements were made. These 
results are consistent with the absence of n loops along 
the grain boundary. 



C. Biepit axial junctions employing CeCh seed layers 



The Ce02 seed layer, as anticipated in section II, may 
produce an artificial GB that can be seen as a result of 
two rotations: a 45° [100] tilt or twist followed by a 45° 
[001] tilt around the c-axis of the (001) film. For this 
junction configuration a d-wave order parameter sym- 
metry would produce 7r-loops, as shown in Fig. 2. We 
notice that such 7r-loops are structurally different from 
those usually obtained by the 45° [001] tilt GB junctions 
based on the traditional biepitaxial and bicrystal tech- 
niques. Due to the microstructure we expect especially 
in the [100] tilt case low critical current densities and 
high normal state resistances. We found that the depo- 
sition conditions to select the uniform growth of YBCO 
45° tilted around the c-axis of the (001) film are critical. 
Preliminary measurements realized on tilt-type junctions 
with a Ce02 seed layer gave evidence of Josephson cou- 
pling in these GBs. The measured IcR-at values are from 
200 fjiV to 750 and are in the typical range of the GBs 
Josephson junctions. 



D. Biepitaxial SQUIDs employing MgO seed layers 



In this section we report on the characterization of dc- 
SQUIDs which are to our knowledge the first employing 
the GBs discussed above. 24 These SQUIDs exhibit very 
good properties, and noise levels which are among the 
lowest ever reported for biepitaxial junctions. 24 Apart 
from implications for applications, these performances 
are important for the study of the transport properties 
of HTS Josephson junctions. In Fig. 7 we show the mag- 
netic field dependence of the voltage at 77 K for different 
values of the bias current for a dc-SQUID with an induc- 
tance of 13 pH. At this temperature IcR-n is about 20 
fiV. The corresponding value of the screening parameter 
/?=2LIc/$o is 0.03. In general low /? values are manda- 
tory to avoid the influence of asymmetric inductances in 
SQUID properties, and this has been crucial for experi- 
ments designed to study the order parameter symmetry. 7 
The presented curves are quite typical. These SQUIDs 
usually work in a wide temperature range from low tem- 
peratures (4,2 K ) up to temperatures above 77 K. The 
maximum working temperature was in this case 82 K. 
The achieved magnetic flux-to- voltage transfer functions 
V$ = dV/d$, where V and $ are the voltage across the 
device and the applied magnetic flux in the SQUID loop 
respectively, are suitable for applications. For instance 
at T = 77 K an experimental value of the SQUID ampli- 
tude voltage modulation AV of 10.4 mV was measured, 
corresponding to V$ = 36.9 //V/<f> 0 - 24 Steps of different 
nature have been recurrently observed in the I-V charac- 
teristics in the washer and hole configurations and char- 
acterized also in terms of the magnetic field dependence 
of the voltage at different values of the bias current. 

The noise spectral densities of the same dc-SQUID 
have been measured at T = 4.2 K and T = 77 K using 
standard flux-locked-loop modulated electronics. The en- 
ergy resolution e = S$/2L (with S$ being the magnetic- 
flux-noise spectral density) at T = 4.2 K and T = 77 K is 
reported in Fig. 8. At T = 4.2 K and 10 kHz, a value of 
S$ = 3 fx$o/y/Hz has been measured, corresponding to 
an energy resolution e = 1.6xl0" 30 J/Hz. This value is 
the lowest reported in the literature for YBCO biepitaxial 
SQUIDs. Moreover, the low frequency 1/f flux noise spec- 
tral density at 1 Hz is more than one order of magnitude 
lower than the one reported for traditional biepitaxials, 
as is also evident from the comparison with data at T — 
4.2 K of Ref. 25 . The lower values of low frequency noise 
are consistent with the absence of 7r- loops on the scale of 
the faceting for these types of GBs, as clearly shown by 
Scanning SQUID Microscopy results. The 7r-loops pro- 
duce some types of spontaneous magnetic flux in the GB 
region, which among other effects tends to degrade the 
SQUID'S noise levels. 14 
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IV. BIEPITAXIAL JUNCTIONS FOR 
EXPERIMENTS ON THE SYMMETRY OF THE 

ORDER PARAMETER AND FOR A 
DEVELOPMENT OF CONCEPTS FOR QUBITS 



The particular junction configurations investigated in 
this work allow some consideration of the possible impact 
of these types of junctions on the study of the Josephson 
effect and the order parameter symmetry in YBCO and 
on the development of concepts for devices. 1 ' 6 ' 9,7 We first 
recall that the biepitaxial technique can provide circuits 
composed completely of junctions without any 7r-loops 
(see Fig. 9a). By varying the interface orientation with 
respect to the [103] electrode orientation, the junction 
properties can be adjusted. On the other hand the tra- 
ditional biepitaxial technique, 11 producing 45° [001] tilt 
GBs (see Fig. 9b) or the types of junctions described in 
the previous section by using Ce02 (see Fig. 9c), can 
controllably generate 7r-loops on macroscopic scales. In 
these schemes we use a corner geometry with a 90° an- 
gle. This angle a can be obviously tuned to enhance 
the effects related to the phase shift (see dashed line in 
Fig. 9b) and this change is particularly easy to realize by 
using the biepitaxial technique. 

In this section we focus our attention mainly on the 
feasibility of the biepitaxial junctions to obtain the dou- 
bly degenerate state required for a qubit. In Ref. ltt the 
design is based on quenching the lowest order coupling 
by arranging a junction with its normal aligned with 
the node of the d-wave order parameter, thus produc- 
ing a double periodic current-phase relation. It has been 
shown that the use of ir phase shifts in a superconducting 
phase qubit provides a naturally bistable device and does 
not require external bias currents and magnetic fields. 16 
The direct consequence is the quietness of the device 
over other designs. A 7r junction provides the required 
doubly degenerate fundamental state, which also mani- 
fests itself in a doubly periodic function of the critical 
current density as a function of the phase. 8 The same 
principle has been used in small inductance five junction 
loop frustrated by a 7r-phase shift. 16 This design pro- 
vides a perfectly degenerate two-level system and offers 
some advantages in terms of fabrication ease and perfor- 
mance. HTS may represent a natural solution for the 
realization of the required 7r-phase shift due to the pair- 
ing symmetry of the order parameter and, therefore, due 
to the possibility of producing ir phase shifts. Experi- 
mental evidence of YBCO 7r -SQUIDs has been given by 
employing the bicrystal technique on special tetracrystal 
substrates. 7 The biepitaxial technique, beyond provid- 
ing junctions with opportune properties, would guaran- 
tee the versatility necessary for the implementation of a 
real device, as shown below. As a matter of fact, we no- 
tice that our technique allows the realization of circuits 
where 7r-loops can be controllably located in part of the 
substrate and separated from the rest of the circuit based 
on "0" -loops, i.e. junctions where no additional 7r phase 



shifts arise. This can be easily made by depositing the 
MgO and Ce02 seed layers on different parts of the sub- 
strate, which will be also partly not covered by any seed 
layer. 

As a test to show how the biepitaxial junctions could 
be considered for preliminary tests and device implemen- 
tation for quantum computing without the topological 
restriction imposed by the bicrystal technique, we refer 
to the structures proposed in Ref^as exemplary circuits. 

The former is composed by a s-wave (S)- d-wave (D)- s- 
wave (S 5 ) double junction connected with a capacitor and 
an ordinary "0" Josephson junction based on s-wave su- 
perconductors (the S-D'-S junction generates the doubly 
degenerate state). The latter consists of a five junction 
loop with a 7r junction. Our technique would combine 
the possibility of placing the ordinary "0" junctions cor- 
responding to the MgO seed layer and to exploit the pos- 
sible doubly degenerate state of asymmetric 45° GB junc- 
tions corresponding to the Ce0 2 seed layer to replace the 
S-D-S' system or the n junction respectively. Our struc- 
ture would be obviously composed only of HTS. In Figs. 
10a and 10b we show how devices for instance such as 
those proposed in Ref. 1 could be obtained by employing 
the biepitaxial technique respectively. The application to 
the five junction loop is straightforward (Fig. 10b) and the 
advantages of this structure have been already discussed 
in Ref. 16 . The biepitaxial technique can offer possible 
alternatives for the realization of the structures above. 
In particular the double junctions of the original S-D- 
S' system can be also replaced by a D'-D-D" structure 
(Fig. 10c) by exploiting our technique, in contrast to the 
bicrystal technology which could not give this possibil- 
ity. Such a configuration could offer some advantages, if 
we consider that asymmetric 45° bicrystal GB Joseph- 
son junctions did not give systematic evidence of the 
doubly degenerate state. The doubly degenerate state 
seems to occur only in high quality low transparency GB 
junctions 8,26 and it is known that S-I-D junctions do not 
have double periodicity of the critical current as a func- 
tion of the phase. 26 A consequence of a possible nodeless 
order parameter 4,23 at the D-D' GB interface could be 
a closer similarity with a S-I-D junction with loss of the 
doubly degenerate state. If this is the case, we spec- 
ulate that the double junctions structure for symmetry 
reasons would produce a leading term in the Josephson 
coupling energy of the form cos20 (double periodic) 
and that the possible dipolar component of the magnetic 
field would be almost completely compensated in this 
configuration. 16 This can be considered as an attempt to 
construct a "microscopic" 20-junction. We finally notice 
that the topological advantages offered by the biepitax- 
ial junctions would therefore be crucial in both the cases 
considered for the realization of the structure in Fig. 10, 
and important to reduce de-coherence effects. Bicrystal 
substrates would in fact impose on the circuit additional 
junctions required by the circuit design and, as a con- 
sequence, generate additional noise and de- coherence in 
the device. 
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V. CONCLUSIONS 

The performance of the presented junctions and 
SQUIDs demonstrates that significant improvements in 
the biepitaxial technique are possible, and the resulting 
devices have potential for applications. We have pre- 
sented a phenomenology that is consistent with the ex- 
pected absence of 7r-loops in 45° [100] tilt and twist grain 
boundaries junctions. The use of a Ce02 rather than a 
MgO seed layer can produce 7r-loops in the same junction 
configurations. The versatility of the biepitaxial tech- 
nique has been recently used to obtain different types 
of grain boundaries. The advantage of placing junctions 
in arbitrary locations on the substrate without imposing 
any restrictions on the geometry, and the ease of obtain- 
ing different device configurations by suitably patterning 
the seed layer, make the biepitaxial technique competi- 
tive for the testing of new concept devices, such as those 
based on 7r-loops. Some simple examples of situations 
in which 7r-loops can be suitably produced in specific 
locations of a more complicated circuit have also been 
discussed. 
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FIG. 1. A schematic representation of the artificial grain 
boundary structure. The boundary is obtained at the inter- 
face between the [001] oriented YBCO film grown on the [110] 
MgO seed layer and the [103] YBCO film grown on the bare 
[110] STO substrate. In contrast with the 45° [001] tilt bicrys- 
tal junctions, in this case the order parameter orientations do 
not produce an additional tt phase shift. 
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FIG. 2. The CeC>2 seed layer produces an artificial GB 
that can be seen as a result of two rotations: a 45° [100] tilt 
or twist followed by a 45° tilt around the c-axis of the (001) 
film. For this junction configuration a d-wave order parameter 
symmetry would produce n -loops. 

FIG. 3. Current vs voltage (I-V) characteristics of the 
biepitaxial junction for temperature close to the critical tem- 
perature. In the inset the I-V curve at T = 4.2 K is shown. 

FIG. 4. a) Scheme of the seed layer patterning, which al- 
lows the measuremment on the same chip of the properties 
of a tilt junction and of junctions whose interface is tilted in 
plane of an angle a = 30°, 45° and 60° with respect the a- or 
b-axis of the (001) YBCO thin film respectively, b) The I-V 
characteristics (measured at T = 4.2 K) of the microbridges 
reported in Fig. 4a. 

FIG. 5. Magnetic-field dependence of the critical current 
of a [100] tilt biepitaxial dc-SQUID. The absolute maximum 
is observed for zero field. A double-period modulation is ob- 
served. The longer period modulation is the diffraction pattern 
due to the magnetic field sensed by a single junction, while 
the shorter period SQUID modulation is shown more clearly 
in the inset (a). In the inset (b) I-V curves are shown as a 
function of an externally applied magnetic field at T = 4.2 K. 
A typical Fraunhofer-like dependence is evident. 



FIG. 6. Scanning SQUID microscope image of a 
200x200^m 2 area along a grain boundary separating a (100) 
region from a (103) region of a thin YBCO biepitaxial film 
grown. The position of the grain boundary is indicated by 
the dashed line. 

FIG. 7. Magnetic field dependence of the voltage of a [100] 
tilt biepitaxial dc-SQUID at 77 K for different values of the 
bias current. 

FIG. 8. Magnetic flux noise spectral densities of a [100] 
tilt biepitaxial SQUID at T=77 K and T=4.2 K. The SQUID, 
with an inductance L=13 pH, was modulated with a standard 
flux-locked-loop electronics. The right axis shows the energy 
resolution. Data at T = 4.2 K are compared with results on 
SQUIDs based on [001] tilt biepitaxial junctions from Ref. 25. 

FIG. 9. a) 3-dimensional view of a SQUID based on 45° 
[100] tilt and twist GBs; no 7T-loops should occur, b) Top 
view of 7T-SQUID based on 45° [001] tilt GBs. c) 3-dimen- 
sional view of a tt-SQUID based on GBs resulting from two 
rotations: a 45° [100] tilt or twist followed by a 45° [001] tilt 

FIG. 10. Scheme of the qubit structure proposed in Ref.l 
designed using the biepitaxial grain boundaries proposed in 
the paper. The double junctions of the original S-D-S' system 
can be also replaced by D'-D-D". 
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We have developed a new way of making grain boundary junctions in YBa 2 Cu 3 0 7 thin films 
by controlling the in-plane epitaxy of the deposited film using seed and buffer layers. 
We produce 45° grain boundaries along photolithographically defined lines. The typical value 
of the critical current density of the junctions is 10 3 -10 4 A/cm 2 at 4.2 K and 10 2 -10 3 A/ 
cm 2 at 77 K, while the rest of the film has a critical current density of 1-3 X 10 6 A/cm 2 at 77 
K. The current-voltage characteristics of the junctions show resistively shunted junction 
behavior and we have used them to fabricate dc superconducting quantum interference devices 
(SQUIDs) which show modulation at temperatures well above 77 K. This is the first 
planar high T c Josephson junction technology that appears readily extendable to high T e 
integrated circuits. 



Most microelectronic applications of the high T c su- 
perconductors will rely upon microbridges such as super- 
conductor-normal metal-superconductor (S-N-S) junctions 
or other weak link geometries instead of superconductor- 
insulator-superconductor (S-I-S) tunnel junctions. A num- 
ber of weak-link structures provide reliable critical current 
reduction in YBa 2 Cu 3 0 7 ( YBCO) thin films; grain bound- 
ary junctions. 1 " 3 traditional edge junctions, 4 microbridges 
grown across a sharp substrate step, 5 and S-N-S-type junc- 
tions with Au, 6 PrBa 2 Cu 3 0 7 , 7 ' 8 or SrTi, ^Nb^Oj (Ref. 9) 
normal layers. Of these various weak-link geometries, grain 
boundary junctions work well at temperatures close to 
T a have reasonably high JcR n products, and show behavior 
explained by a resistively shunted junction model with a 
relatively uniform current density. To date, the high angle 
grain boundary junctions that have been used were formed 
either by fusing differently oriented SrTi0 3 substrates 1 or 
by patterning the randomly occurring grain boundaries in 
granular films. 2 ' 3 Superconducting quantum interference 
devices (SQUIDs) made from the grain boundaries on Sr- 
Ti0 3 bicrystals had reasonable signal and low noise 10 up to 
temperatures very close to the 7V of the film. However, this 
SrTi0 3 bicrystal technique has a major drawback: it cannot 
be readily extended to integrated circuits. On the other 
hand, SQUIDs made from granular films are plagued by 
low yield and the presence of grain boundaries in the 
SQUID loop itself. These lead to excessive flux noise and 
hysteresis in the voltage-flux response of the SQUID. 

In this letter, we report for the first time the successful 
fabrication of high T c weak-link Josephson junctions oper- 
ating at and above 77 K that are made from grain bound- 
ary junctions whose locations are determined by conven- 
tional photolithography. This result demonstrates the 
feasibility of an extendable technology for fabricating high 
T c Josephson devices which avoids the problems associated 
with the devices made from granular thin films and on 
bicrystal substrates. We reproducibly fabricate 45° grain 
boundaries in YBCO films by controlling their in-plane 
epitaxy using seed and buffer layers deposited on r-plane 
sapphire substrates, a method we call "bi-epitaxy." The 
method is quite general, and can be used on a wide variety 
of substrates* 



Recently, we reported 11 the existence of two in-plane 
epitaxial orientations of c-axis oriented epitaxial YBCO 
thin films on yttria-stabilized zirconia (YSZ), resulting 
from the poor lattice match between YBCO and YSZ. This 
observation led us to realize that it would be possible to 
make bi-epitaxial grain boundaries on such a substrate us- 
ing seed and buffer layers, or, in an alternative geometry, 
on a well lattice-matched substrate using poorly lattice- 
matched seed layers. For example, we previously re- 
ported 12 the growth of YBCO thin films on r-plane sap- 
phire using buffer layers such as SrTi0 3 , CaTi0 3 , YSZ, and 
MgO. We emphasized the importance of the in-plane epi- 
taxy of YBCO, particularly because high-angle grain 
boundaries in YBCO films behave as weak links, leading to 
low critical current densities and high surface resis- 
tances. 12 * 13 Our efforts towards minimizing the occurrence 
of these high-angle grain boundaries taught us how to con- 
trol them and led us to design and demonstrate a variety of 
45° weak-link junctions. One example takes advantage of 
the differing epitaxial directions of two different layers, 
MgO and SrTi0 3 , when they are deposited on an r-plane 
sapphire substrate. Using these two layers as a seed and a 
buffer layer, we have succeeded in fabricating bi-epitaxial 
45° grain boundaries in YBCO on sapphire substrates by 
the following process. 

Schematic views of our device structure are shown in 
Fig. 1. We first deposit 3-30 nm of epitaxial MgO as a seed 
layer on an r-plane sapphire substrate. MgO is produced by 
laser ablating a Mg metallic target in a 2-20 mTorr oxygen 
atmosphere with a substrate temperature of 300-700 °C. 
We then mask the MgO with conventional photoresist and 
remove it from, for example, half of the substrate by either 
Ar ion beam milling or chemical wet etching with dilute 
phosphoric acid. We then grow 10-100 nm of epitaxial 
SrTi0 3 buffer layer by laser ablation on both the exposed 
sapphire surface and on the patterned MgO seed layer. The 
deposition conditions are an oxygen pressure of 100-200 
mTorr and a substrate temperature of 710-760 °C The 
SrTi0 3 film grows in two different orientations separated 
by a 45° grain boundary. In the growth plane the epitaxial 
relations are SrTi0 3 [l 10]//Al 2 O 3 [l 120] and 
SrTiO 3 [100]//MgO[100J//Al 2 O 3 [1120], as illustrated in 
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FIG. 1. Schematic view of the device structure with an MgO seed layer 
and a SrTiOj buffer layer on an r-plane sapphire substrate, (a) top view, 
(b) side view. 



Fig. 1. We then immediately deposit YBCO, which grows 
epitaxially everywhere and thereby reproduces the grain 
boundary in the SrTi0 3 buffer layer. Finally the YBCO is 
patterned into an appropriate geometry. 

The virtue of this structure is that, apart from the 
predefined 45° grain boundaries, there are no other high- 
angle grain boundaries in the YBCO film. It has been re- 
ported 12 that YBCO films grown on SrTiOj buffer layers 
on /--plane sapphire substrates have high critical current 
density and low surface resistance. In addition, the SrTiOj 
layer grown on the MgO buffer layer does not have any 
high-angle grain boundaries, as demonstrated in a <f> scan of 
the SrTi0 3 (101), although it is found 13 that YBCO de- 
posited on MgO may have a number of high-angle grain 
boundaries. The probable explanation for this effect is that 
the energy needed to nucleate 45° misoriented SrTi0 3 
grains on MgO is much higher than in the case of YBCO 
grown on MgO. Hence the YBCO film on SrTi0 3 on MgO 
also lacks any high-angle grain boundaries and has a high 
critical current density. Microbridges formed either side of 
the predetermined grain boundary have critical current 
densities of 1-3 X 10 6 A/cm 2 at 77 K. In order to show that 
the critical current reduction is indeed due to the grain 
boundary rather than simply the step at the edge, we 
etched only part of the MgO, reducing its thickness on one 
side from 20 to 10 nm, such that YBCO/SrTi0 3 /MgO 
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FIG. 2. X-ray <b scan of YBCO ( 103) of the unpatterned structure. The 
peaks labeled with "A" belong to YBCO/SrTKVAUO, and the peaks 
labeled with "B" come from YBCO/SrTiOyMgO/AUO^. 



remained on the sapphire substrate on both sides of the 
step.The critical density across this step remained larger 
than 10 6 A/cm 2 at 77 K. 

A 0 scan of the YBCO ( 103) before the last patterning 
process is shown in Fig. 2. The peaks occurring 45° apart 
indicate that half of the YBCO is rotated 45° relative to the 
other half. Suitable etching experiments have shown that 
the peaks labeled with "A" belong to YBCO/SrTi0 3 / 
Al 2 Oj and the peaks labeled with "B" come from YBCO/ 
SrTi03/MgO/Al 2 0 3 . Cross-sectional transmission electron 
microscope images and microdiffraction patterns show that 
the 45* grain boundary begins at the end of the MgO seed 
layer. 14 

As the first example of a device using multiple junc- 
tions made with our process we fabricated a square washer 
SQUID 15 as shown in the inset of Fig. 3. It was patterned 
across the grain boundary by photolithography and wet 
etched with dilute phosphoric acid. The current-voltage 
characteristic (7-FC) of the SQUID at 4.2 K, which ex- 
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FIG. 3. Current (vertical) -voltage (horizontal) characteristics of a dc 
SQUID at 4.2 K. The inset shows the geometry of the SQUID. 
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FIG. 4. Modulation voltage of the current-biased dc SQUID vs applied 
magnetic flux at various temperatures. 

hibits resistively shunted Josephson junction behavior, is 
shown in Fig. 3. Each junction was 4 [im wide and 0.2 /xm 
thick. The critical current of about 10 jxA at 4.2 K trans- 
lates to about 10 3 A/cm 2 critical current density across the 
junctions. The 7^ rt product of the junctions is about 420 
fiV at 4.2 K. The l-VC of the SQUID at 77 K shows no 
zero resistance part, although it is nonlinear. The small 
junction coupling energy I/^q/Itt compared to the thermal 
energy k B T may lead to thermally activated phase slippage 
across the junctions, resulting in resistance at all cur- 
rents. 16 Detailed transport properties of the junctions in- 
cluding the magnetic field dependence of the critical cur- 
rent will be reported elsewhere. 17 

The voltage across the current-biased dc SQUID as a 
function of applied field is shown for various temperatures 
in Fig. 4. The data were taken in a bandwidth of 0-10 Hz 
without any signal averaging. The observed period corre- 
sponds to the expected value from the geometry of the 
SQUID. We believe that the reduction of the modulation 
voltage at high temperatures is due to the decreasing crit- 
ical current and the large inductance of the SQUID. We 
have observed dc SQUID operation at temperatures as 
high as 88 K and have operated such a SQUID in a flux- 
locked loop up to 83 K. The noise of this SQUID is similar 
to that reported in dc SQUIDs made using bicrystal sub- 
strates. 10 Detailed performance of the SQUID will be pub- 
lished elsewhere. 18 

A further example of small-scale integration using 
biepitaxial grain boundary junctions is provided by a flux 



shuttle 19 we have fabricated. It involves two dc SQUIDs as 
sensors and eleven more bi-epitaxial junctions as switching 
elements. Our work on such integrated circuits is continu- 
ing and will be reported later. 

In summary, we have reported the development of bi- 
epitaxial grain boundary junctions in YBa 2 Cu 3 0 7 at mul- 
tiple and predetermined locations by using only standard 
photolithographic techniques to control in-plane epitaxy 
with seed and bufTer layers. The junctions have current- 
voltage characteristics that are well described by the resis- 
tively shunted junction model and dc SQUIDs fabricated 
using these junctions show modulation at temperatures as 
high as 88 K. By increasing the critical current density of 
the junctions and lowering the inductance of the SQUID, 
the performance of the SQUIDs should improve further, 
especially at 77 K. We believe that this technology, based 
on bi-epitaxial grain boundaries, can be extended to inte- 
grated circuits in the near further, and we are planning to 
demonstrate such circuits. 

We would like to thank John Rowell, Mac Beasley, 
Ted Geballe, Bob Hammond, John Clarke, Aharon Kapit- 
ulnik, Randy Simon, and Roger Barton for helpful discus- 
sions on this work. We also would like to thank JefF Rosner 
at Hewlett-Packard for the TEM studies. 
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Anomalous Periodicity of the Current-Phase Relationship of Grain-Boundary 
Josephson Junctions in High-T c Superconductors 
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The current-phase relation (CPR) for asymmetric 45° Josephson junctions between two d-wave 
superconductors has been predicted to exhibit an anomalous periodicity. We have used the single- 
junction interferometer to investigate the CPR for this kind of junctions in YBa 2 Cu 3 0 7 -x thin films. 
Half-fluxon periodicity has been experimentally found, providing a novel source of evidence for the 
d-wave symmetry of the pairing state of the cuprates. 
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There is growing evidence in favor of the d^^-wave 
symmetry of the pairing state of the high-temperature 
superconductors. 1 An unconventional pairing state re- 
quires the existence of zeros of the order parameter in 
certain directions in momentum space. Thermodynamic 
and spectroscopic measurements do indeed suggest their 
existence, but by themselves they do not exclude conven- 
tional 5- wave pairing with nodes. 1 Direct evidence for 
the d-wave pairing state is provided by phase-sensitive 
experiments, which are based on the Josephson effect. 2 
Quite generally, the current-phase relationship (CPR) of 
a Josephson junction, J(y>), is an odd periodic function 
of <p with a period 2tt. 3 Therefore I(ip) can be expanded 
in a Fourier series 



I (if) = h sin <p + h sin 2(p + 



(i) 



In the tunnel limit we can restrict ourselves to the first 
two terms in Eq. (1). Since the order parameter is 
bound to the crystal lattice, I (if) of a weak link de- 
pends on the orientation of the d-wave electrodes with 
respect to their boundary. The existing phase-sensitive 
experiments exploit possible sign changes of h between 
different geometries. 2 In this Letter we present a new 
phase-sensitive experimental test of the symmetry of the 
pairing state of the cuprates. Namely, in certain geome- 
tries, the ii term should vanish by symmetry. In such 
cases, the CPR should exhibit an anomalous periodicity. 

Let us analyze the angular dependence of hi a junc- 
tion between two macroscopically tetragonal d-wave su- 
perconductors. As emphasized in Ref. 4, also heavily 
twinned orthorhombic materials such as YBa2Cu307_ x 
belong to this class, if the twin boundaries have odd sym- 
metry. We consider first an ideally flat interface between 
the superconducting electrodes. Let 9\ (0 2 ) denote the 
angle between the normal to the grain boundary and the 
a axis in the electrode 1 (2), see Fig. 1. If we keep only 
the lowest-order angular harmonics, the symmetry of the 
problem dictates that 4 

h = h cos 20i cos 20 2 + Is sin 29 x sin 20 2 . (2) 



The coefficients I C ,I 3 are functions of the barrier 
strength, temperature T, etc. The I 2 term results from 
higher-order tunneling processes and we neglect its weak 
angular dependence. It is seen from Eq. (2) that the cri- 
terion for the observation of an anomalous CPR, h = 0, 
is realized for an asymmetric 45° junction, i.e. a junction 
with 6 1 — 45° and 0 2 = 0. For an interface which is not 
ideally flat, 0* = 0i(x) are functions of the coordinate x 
along the junction, h = 0 remains valid even in this 
case, if the average values ($i(x)) = 45° and (0 2 (x)) = 0. 
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FIG. 1. Sketch of the interferometer (not in scale). 

The I 2 term is present also in weak links based on 
conventional s-wave superconductors but for all known 
types of weak links |J 2 /Ji| < 1. For instance, for a tunnel 
junction \I 2 /h\ 1. For a SNS junction, / oc sinv?/2 
at T = 0, 5 and the Fourier expansion Eq. (1) leads to 
hlh = -2/5. Therefore, a possible experimental ob- 
servation of \I 2 /h\ > 1 in an asymmetric 45° junction 
provides direct evidence of d-wave symmetry of the pair- 
ing state in the cuprates. 

We have investigated the CPR of YBa 2 Cu 3 0 7 -x 
thin film bicrystals with asymmetric 45° [001]-tilt 
grain boundaries 6 " 8 as sketched in Fig. 1, using a 
single-junction interferometer configuration in which the 
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Josephson junction is inserted into a superconducting 
loop with a small inductance L. In a stationary state 
without fluctuations, the phase difference <p across the 
junction is controlled by applying external magnetic flux 
$ e penetrating the loop: tp = ip e - (3f{ip). Here <p e = 
27r<I> e /$o is the external flux normalized to the flux quan- 
tum $ 0 = 2.07 x 10 -15 Tm 2 . The CPR is written as 
I (if) = Iof{(f), where h is the maximal Josephson cur- 
rent. P = 2ttLIo/^o is the normalized critical current. 
In order to obtain the CPR for the complete phase range 
— 7r < if < 7r the condition (3 < 1 has to be fulfilled, 
because for ft > 1 the curve <p((p e ) becomes multivalued 
and there are jumps of tp and a hysteresis for a sweep 
of (f e - Following Ref. 3, we express the effective induc- 
tance of the interferometer using the derivative /' with 
respect to tp as L int = £[1+1 //'(¥>)]■ The inductance can 
be probed by coupling the interferometer to a tank cir- 
cuit with inductance L>t, quality factor Q, and resonance 
frequency t^o- 9 External flux in the interferometer is pro- 
duced by a current Idc + Irf i n the tank coil and can be 
expressed as v? e = 27r(/d C -h/ r /)M/<I>o = Vdc + Vrf, where 
M 2 = k 2 LL,T, and k is a coupling coefficient. Taking into 
account the quasiparticle current in the presence of volt- 
age V across the junction the phase difference is given 
by the relation tp = tpdc + <Prf - 0f(<p) ~ 27rr(ip)V/^o, 
where r(ip) = L/Rj(ip) and Rj{<p) is the resistance of 
the junction . In the small-signal limit <p r f < 1 and in 
the adiabatic case ujt <IC 1, keeping only the first-order 
terms, the effective inductance L e ff of the tank curcuit- 
interferometer system reads 



*>eff 



Thus the phase angle a between the driving current and 
the tank voltage U at the resonant frequency of the tank 
circuit ujq is 



tana(<p) = 



k 2 Q/3f'(<p) 
l+0f'(f)' 



(3) 



Using the relation [1 +Pf'{ip))dy> = d<pd c valid for <p r f < 
1 and wr < 1, one can find the CPR from Eq. (3) by 
numerical integration. 

The advantage of the measurement of the CPR of an 
asymmetric 45° junction with respect to the by- now stan- 
dard phase-sensitive tests of pairing symmetry based on 
the angular dependence of I\ is twofold. First, it avoids 
the complications of the analysis of experiments caused 
by the presence of the term 7 S . 4 Second, a flux trapped 
in the SQUID does not invalidate the conclusions about 
the ratio \h/h\ and hence about the pairing symmetry, 
while this is not the case in standard phase-sensitive tests 
of the rf-wave symmetry of the pairing state. 10 

The films of thickness 100 nm were fabricated us- 
ing standard pulsed laser deposition on (001) oriented 
SrTiC>3 bicrystalline substrates with asymmetric [001] tilt 
misorientation angles 45° ± 1°. They were subsequently 
patterned by Ar ion-beam etching into 4x4 mm 2 square 



washer single-junction interferometer structures (Fig. 1). 
The widths of the junctions were 1 -=- 2 fim. The washer 
square holes had a side-length of 50 /zm. This geometry 
of the interferometer gives L w 80 pH. The resistance of 
the junction is higher than 1 O and the condition for the 
adiabatic limit ljt <fC 1 is satisfied. For measurements of 
a(<pdc), several tank circuits with inductances 0.2 -=- 0.8 
/iH and resonance frequencies 16 4- 35 MHz have been 
used. The unloaded quality factor of the tank circuits 
70 < Q < 150 has been measured at various tempera- 
tures. The coupling factor k was determined from the 
period Aide of <x{hc) using MAIdc — $o- Its value var- 
ied between 0.03 and 0.09. The amplitude of I r f was 
set to produce the flux in the interferometer lower than 
0.1*0- 

The measurements have been performed in a gas- flow 
cryostat with a five-layer magnetic shielding in the tem- 
perature range 4.2 < T < 90 K. The experimental setup 
was calibrated by measuring interferometers of the same 
size with 24° and 36° grain boundaries. We have stud- 
ied 5 samples, out of which sample No. 1 exhibited the 
most anomalous behavior. Samples Nos. 2,3 were less 
anomalous and the remaining two samples had high crit- 
ical currents and their I((p) was conventional. In Fig. 2 
we plot the phase angle a as a function of the dc current 
Idc for samples Nos. 1,2. The data for the 36° junction 
is shown for comparison. Note that at T = 40 K the 
periodicity of a(Id C ) is the same for all samples. 
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FIG. 2. Left panel: Phase angle between the driving current 
and the output voltage measured for the sample No. 1 at dif- 
ferent temperatures as a function of the dc current Idc The 
curves are shifted along the y axis and the data for T = 30 
and 40 K are multiplied by factor 4 for clarity. From top to 
bottom, the data correspond to T =4.2, 10, 15, 20, 30 and 40 
K. The data measured on 36° bicrystals (0i « 36° , 62 « 0) 
at T =40 K in the same washer geometry are shown for com- 
parison (open circles). Right panel: The same for the sample 
No. 2. From top to bottom, the data correspond to T =4.2, 
10, 15 and 20 K. 

We assume that the period of a(Id c ) at T =40 K, 
A/dc =9.6 £tA, corresponds to Aipdc — 27r. In order 
to determine the CPR we take <pd c = 0 at a maximum 
or minimum of a. This is necessary in order to satisfy 
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J(ip — 0) = 0, as required by general principles. 3 The ex- 
perimentally observed shift of the first local extreme of 
a{I dc ) from he = 0 (Fig. 2) can be due to flux trapped 
in the interferometer washer. Most probably, this flux re- 
sides in the long junction of the interferometer. The long 
junction does not play an active role because the Joseph- 
son penetration depth is much shorter than its length, 
and external fields producing I dc are smaller than its crit- 
ical field. Nevertheless the long junction sets the phase 
difference for he = 0 at the small junction. 
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FIG. 3. a) Josephson current through the junction for the 
sample No. 1 as a function of the phase difference y>, deter- 
mined from the data in Fig. 2. The statistics of a{(p) was 
improved by folding the data back to the interval {0, 7r) and 
taking an average. The symmetry a(tp) — a(-ip) was as- 
sumed, b) The same for the sample No. 2. The data for the 
asymmetric 36° bicrystal at T = 40 K (open circles) is also 
shown. 

In Fig. 3, we show the CPR determined from the data 
in Fig. 2. For all curves we have performed a minimal 
necessary shift consistent with I [op = 0) = 0. Thus, for 
the samples Nos. 1,2 we have assumed that at ipdc — 0 
a minimum of a(<pd c ) is realized. For an interferometer 
with a conventional s-wave weak link (and also for the 36° 
junction), at ifdc = 0 a maximum of a((fdc) is realized. 
Note that the maximum (minimum) of a(<pd c ) at (fd c = 0 
implies a diamagnetic (paramagnetic) response of the in- 
terferometer in the limit of small applied fields. In Fig. 4 
we plot the coefficients h and h determined by Fourier 
analysis of the CPR for the sample No. 1 at various tem- 
peratures. With decreasing T, |/ 2 | grows monotonically 
down to T — 4.2 K, while the I\ component exhibits only 
a weak temperature dependence. 



Our experimental results can be understood as follows. 
Deviations from ideal geometry of the asymmetric 45° 
junction, (0i> = 45° + «i and (9 2 ) = <*2, lead to a fi- 
nite value of h. Thus, imperfections of the junction 
increase its critical current. For this reason we believe 
that samples Nos. 2-5 contain imperfections and from 
now we concentrate on nearly ideal junctions (such as 
sample No. 1) with |ai|,|a 2 | < 1. For such junctions, 
the ratio I 2 /h exhibits the following temperature de- 
pendence. For T — 0, \hlh \ » 1. The region T ~ T c 
can be analyzed quite generally within Ginzburg-Landau 
theory. Let the electrodes be described by (macro- 
scopic) order parameters Ai, 2 = |A|e iv31 ' 2 . Then the 
phase-dependent part of the energy of the junction is 
E = a[Ai + H.C.] + 6[(Ai A£) 2 + H.C.] + . . . where 
a, 6, . . . depend weakly on T. 11 Thus for T close to T c we 
estimate h a |A| 2 oc (T c -T) and J 2 a |A| 4 <x (T c -T) 2 , 
leading to \h/h\ < 1- These expectations are qualita- 
tively consistent with the experimental data shown in 
Fig. 4. 
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FIG. 4. Temperature dependence of the Fourier expansion 
coefficients Ii,2 determined from the experimental data in 
Fig. 3a. Solid lines are the Fourier expansion coefficients for 
the numerical data in Fig. 5. 

So far, our discussion was based solely on symmetry 
arguments. Let us attempt a more quantitative analy- 
sis of our data now. Two different microscopic pictures 
of asymmetric 45° Josephson junctions between d-wave 
superconductors have been considered in the literature. 
The first picture assumes a microscopically tetragonal 
material and an ideally flat interface. 11 " 13 Within this 
picture, there are two contributions to the Josephson 
current. The first is due to bulk states and in the tun- 
nel limit it is well described by the Sigrist-Rice term 
I c in Eq. (2). 14 The second is due to mid-gap states 
which develop close to the surfaces of unconventional 
superconductors. 15 I{<p) for the sample No. 1 calculated 
according to the model of Ref. 12 is shown in Fig. 5. 
The experimental data can be fitted by a relatively broad 
range of barrier heights. However, if we require the 36° 
junction to be fitted by the same (or smaller) barrier 
height as for the 45° junction, we conclude the barrier 
must be rather low. 16 The T dependence of I(<p) requires 
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a choice of T c « 60 K in the non-selfconsistent theory 
of Ref. 12. The reduction from the tyulk T c = 90 K is 
probably due to a combined effect of surface degradation 
and order-parameter suppression at the sample surface. 
The temperature dependence of the ratio of the 7r and 
27r periodic components in I(tp) is seen to be in qualita- 
tive agreement with experimental data in Fig. 3a. This 
is explicitly demonstrated in Fig. 4 where we compare 
the experimentally obtained 7i, 2 with the results of the 
Fourier analysis of the curves in Fig. 5. The divergence of 
7 2 as T — ► 0 is an artifact of the ideal junction geometry 
assumed in Ref. 12. If the finite roughness of the inter- 
face is taken into account, this divergence is cut off and 
the experimental data in Fig. 4 do indeed resemble the- 
oretical predictions for a rough interface. 13 However, the 
nonselfconsistent theory of Ref. 12 is unable to explain 
the experimentally observed steep CPR close to the min- 
ima of the junction energy (see Fig. 3a). In the limit of 
vanishing barrier height, the theoretical CPR does have 
steep portions, but these are located close to the maxima 
of the junction energy (see also Ref. 11). 
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FIG. 5. I(<p) calculated according to Eq. (64) of Ref. 12 for 
a junction with 0 X = 45.5° , 6 2 = 0, Xd = 1.5, k = 0.5, and 
T c = 60 K. I(<p) at T = 40 K for the 36° bicrystal (open 
circles) was calculated for the same parameters except for 
0i = 36°. 

In a different approach to the asymmetric 45° junction, 
one assumes a heavily twinned orthorhombic material 
(which is macroscopically tetragonal, however) and/or 
a meandering interface with 0* - 0i(x). 17)18 Hence the 
critical current density j c (x) is a random function with 
a typical amplitude (|j c (z)l) ~ 3c- If the average criti- 
cal current along the junction (j c ) < jc, a spontaneous 
flux is generated in the junction, and I/2//1I » 1. 17,18 In 
particular, for = 45° and (0 2 ) = 0, there is an equal 
amount of parts having positive and negative j c , leading 
to (Jc) = 0 and h = 0. Note that also within the picture 
of Refs. 17,18, the d-wave symmetry of the pairing state 
is crucial, otherwise the condition (j c ) <§; j c is difficult to 
satisfy. 



Our present understanding of I{<p) in the asymmetric 
45° junction is only qualitative. There is considerable 
experimental evidence 6 " 8 that the grain boundary junc- 
tions are at most piecewise flat. However, we cannot say 
whether the shape of I{<p) is dominated by the mid-gap 
states in the microscopically flat regions, or by sponta- 
neous flux generation due to the spatial inhomogeneity 
of the junction- This issue requires further study. 

In conclusion, we have measured the magnetic field re- 
sponse of a single-junction interferometer based on asym- 
metric 45° grain-boundary junctions in YBa2Cu307_a; 
thin films. Half-fluxon periodicity has been experimen- 
tally found, in agreement with theoretical predictions for 
d a .2_ y 2-wave superconductors. Hence, our results provide 
a novel source of evidence for the d-wave symmetry of the 
pairing state in the cuprates. 
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10. I declare that, based on my experience in the field of bi-epitaxial technology, 
neither the bi-epitaxial technology described in Char et al. nor the best quality crystal 
structures available for bi-epitaxial Josephson junction technology were sufficiently advanced 
at the time of filing of the Application to prepare a clean Josephson junction such as that 



described in the Application. 

11. I further declare, under penalty of perjury under the laws of the United States 
of America, that all statements made herein of my own knowledge are true and that these 
statements were made with the knowledge that willful false statements and the like are 



punishable by fine or imprisonment, or both, under Section 1001 or Title 18 of the United 
States Code. 





